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Abstract 

Standard perturbative (or mean field theory) techniques are not adequate for 
studying the finite-temperature electroweak phase transition in some cases of 
interest to scenarios for electroweak baryogenesis. We instead study the prop- 

(N 

\ erties of this transition using the renormalization group and the e-expansion. 

\ This expansion, based on dimensional continuation from 3 to 4— e spatial di- 

\ mensions, provides a systematic approximation for computing the effects of 

I ' (near)-critical fluctuations. The e-expansion is known to predict a first-order 

^ \ transition in Higgs theories, even for heavy Higgs boson masses. The validity 
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of this conclusion in the standard model is examined in detail. A variety of 



5h ' physical quantities are computed at leading and next-to- leading order in e. For 

moderately light Higgs masses (below 100 GeV), the e-expansion suggests that 
the transition is more strongly first order than is predicted by the conventional 
analysis based on the one-loop (ring-improved) effective potential. Neverthe- 
less, the rate of baryon non-conservation after the transition is found to be 
larger than that given by the one-loop effective potential calculation. Detailed 
next-to-leading order calculations of some sample quantities suggests that the 
e-expansion is reasonably well behaved for Higgs masses below 100-200 GeV. 
We also compare the e-expansion with large- results (where N is the number 
of scalar fields) and find that the e-expansion is less well behaved in this limit. 
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I. INTRODUCTION 



The electroweak phase transition has received considerable attention in the past few 
years because of its role in electroweak scenarios for baryogenesis [Q. Important problems 
include the determination of the order and strength of the phase transition as a function of 
the Higgs boson mass, and the calculation of the rate of baryon non-conservation near the 
transition. In most previous work, the principle tool for studying these problems has been 
the one-loop, ring-improved, finite-temperature effective potential for the Higgs field.Q As 
discussed below, the use of the ring-improved loop expansion appears untrustworthy in many 
applications, including the Higgs mass bounds for baryogenesis in the minimal standard 
model The failure of the loop expansion is exactly analogous to the failure of mean 

field theory for describing critical behavior in condensed matter systems. One alternative 
method for studying critical systems, which has had considerable success in a variety of 
theories, is to use the e-expansion to improve the organization of the perturbation series. 
The e-expansion generalizes the three spatial dimensions of the theory to 4— e dimensions, 
solves the theory when e is small, and then extrapolates to e = 1 0]. In this paper, we 
shall apply the e-expansion to the electroweak phase transition with the goal of estimating 
various parameters of the phase transition. We compute a number of observables at or near 
the transition including the scalar correlation length, latent heat, surface tension, free energy 
difference, bubble nucleation rate, and the sphaleron transition (or baryon non-conservation) 
rate. Next-to-leading order results are presented for the correlation length and the latent 
heat. We discuss in detail the expected reliability of the e-expansion results by (i) comparing 
leading and next-to-leading order contributions, (ii) comparing results of the e-expansion 
with standard perturbation theory in the limit of a small Higgs mass (where the ordinary 
loop expansion remains reliable), and (iii) comparing e-expansion results when the number 



^ "Ring-improved" refers to the inclusion of lowest-order Debye screening corrections in the bare gauge and 
Higgs field propagators. 
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of scalar fields is large to direct large- A?" calculations. 

In the remainder of the introduction, we review the need for detailed knowledge of the 
electroweak phase transition in order to determine the viability of electroweak baryogenesis 
scenarios, and discuss the reliability of the standard loop expansion. Then we briefiy sketch 
the rather successful e-expansion results for critical behavior in simple 4>'^ theories, present an 
overview of the application of the e-expansion to the standard electroweak theory, and discuss 
some possible pitfalls of this approach. We also review the existence of a non-perturbative 
magnetic mass in SU(2) theory and discuss its implications for the e-expansion. Detailed 
calculations follow in subsequent sections, which will be outlined later in this introduction. 

A. Motivation 

One of Sakharov's three famous requirements for baryogenesis is that the universe be 
out of equilibrium. In electroweak scenarios, baryogenesis occurs during the electroweak 
phase transition and Sakharov's condition is met if that transition is first order rather than 
second order. First-order transitions proceed by the nucleation, expansion, and coalescence 
of bubbles of the new phase inside the old — a highly non-equilibrium process. However, the 
phase transition must not merely be first order, it must be sufficiently strongly first order. 
In particular, one constraint arises from considering the rate of baryon number violation at 
the end of the phase transition, when baryogenesis is completed. This rate is exponentially 
sensitive to the expectation value {(p) of the Higgs field responsible for symmetry breaking, 
behaving as Fb ~ exp{ — (const. where /3 is the inverse temperature and g the elec- 
troweak coupling. If {4>) is too small after the phase transition, then baryon-number violation 
continues unabated after baryogenesis and the universe relaxes back to an equilibrium state 
of zero net baryon number. Any successful scenario requires that {(f>) be sufficiently large 
at the end of the transition so that baryon number violation is effectively shut ofi:. Thus, 
one important goal of studies of the electroweak phase transition has been to determine the 
strength of the transition, as measured by the size of the discontinuity in (0) across the 
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transition. 

These issues have been previously studied using the one-loop ring-improved finite- 
temperature potential for the Higgs field With this approximation to the potential, 
one finds a first-order phase transition and the Higgs expectation value immediately after 
the transition can easily be evaluated. In particular, one may study the strength of the tran- 
sition in the simplest interesting model: the minimal standard model with a single Higgs 
doublet. This is presumably only a toy model because, unlike multiple-Higgs models, it is 
unclear whether the minimal model incorporates sufficient CP violation for baryogenesis.0 
Nevertheless, because of its small set of unknown parameters, the minimal model provides a 
relatively uncluttered testing ground for techniques used to analyze the phase transition. In 
this model, employing the one-loop ring-improved potential. Dine et al. found that baryon 
number violation is turned off at the end of baryogenesis only if the zero-temperature Higgs 
mass mh(0) is less than 35 to 40 GeV. A successful scenario for electroweak baryogenesis is 
then incompatible with the experimental limit mh(0) > 60 GeV. 

This disaster is avoidable in multiple-Higgs models, but in any model it will be impor- 
tant to understand the reliability of predictions based on the one-loop potential. So let 
us stick with the minimal model and review the reliability of the above constraints. The 
loop expansion parameter at high temperature is not simply but rather g'^T /M, where 
the characteristic mass scale is set by the ly-boson mass (at temperature T), M gcj). 
Computing M at the phase transition, one finds that the loop expansion parameter g^T jM 
is of order Xj g^ . (A review of this power counting may be found in section 2 of Ref. 0.) 
Hence, the loop expansion breaks down, and calculations based on the one-loop potential are 
unreliable, unless \ <^ g^] that is, unless the Higgs mass mh(0) at zero temperature is much 



^There have been several attempts to study corrections to these resuhs using "super-daisy" or similar 
techniques based on solving mass- gap equations with partial inclusion of higher order terms. For a criticism 
of such methods, see Appendix]^ of Ref. ||^ and the discussion of overlapping momenta in Ref. ||^. 

^For the optimistic assessment, see Ref. [^. 
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less than the W-mass Mw(O).0 This criteria is not obviously satisfied by the mh(0) bound 
of 35 to 40 GeV found by Dine et al, and the importance of higher-order corrections will be 
determined by all the detailed numerical factors left out of this facile estimate of the loop 
expansion parameter. To test the reliability of the one-loop results, two-loop corrections to 
the potential were computed in Refs. and are displayed at the critical temperature in 
Fig. ID for mh(0) = 35 GeV and a top quark mass mt(0) = 110 GeV. Computed in Landau 
gauge, (0) increases by only 20% when two-loop corrections are included. But the height of 
the barrier in the potential increases by a factor of almost three! This dramatic change in 
the height suggests that the loop expansion is unreliable.0 One might optimistically hope 
that the correction to the logarithm of the baryon violation rate Fb will be small because 
In Fb is proportional to (0) , which has a relatively small two-loop correction. However, this 
proportionality is only valid at leading order, and IuFb has its own, independent correc- 
tions controlled by X/g"^. There is no obvious reason why these corrections might not be 
numerically large, similar to the correction to the hump in the potential. 




(|)[GeV] 

FIG. 1. The effective potential at the critical temperature for TOh(O) = 35 GeV and TOt(O) = 110 GeV. 
The dashed and solid lines are the one-loop and two-loop results respectively. 

Because of these problems with the loop expansion, it is clearly appropriate to explore 

''With our conventions M^(0) = .9^((/')V4 and ml{Q) = A (0)2/3, so that A/g^ = 3ml(Q)/UI^{0). 

^Ref. 1^ missed this signal of unreliability because, in numerical work, only the correction to {(/)) was 
computed. 
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other techniques, such as the e-expansion, for studying the electroweak phase transition. 



B. A success story: pure scalar theory. 



Before discussing the electroweak theory, it will be useful to review briefly the e-expansion 
in a much simpler model where it enjoys some success. Consider a real scalar field theory 
in 3 + 1 spacetime dimensions at temperature T, with the Euclidean Lagrangian 



One could attempt to analyze the phase transition by computing the one-loop ring-improved 
potential. With this approximation one would predict a first-order phase transition. How- 
ever, this computation also yields a loop expansion parameter XT/mcs{T) of order 1 at the 
transition, and so the one- loop result cannot be trusted.^ In fact, this model is in the same 
universality class as the Ising model and is known to have a seconc^-order phase transition. 

The e-expansion analysis of this transition is well known and proceeds as follows. The 
study of a second-order phase transition requires exploring the infrared behavior of the 
theory. In the Euclidean formulation of equilibrium finite-temperature quantum field theory, 
space-time is periodic in the time direction with period (3 = 1/T. For spatial distances r large 
compared to jS, fluctuations in the Euclidean time direction are irrelevant and the theory 
reduces to an effective three-dimensional theory. Up to corrections suppressed by powers of 
/3/r, the only effect of decoupled temporal fluctuations is to renormalize the parameters of 
the three dimensional theory, i.e., masses and coupling constants, from their original four- 
dimensional values. The relationship between the three- and four-dimensional parameters 
is perturbatively calculable. In the case at hand, one finds that the mass parameters in the 
three- and four-dimensional theories are related by 






^Again, see section 2 of ref. [|j for a brief review in the same language used here. 
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The second term is the standard one-loop thermal contribution to the mass and is responsible 
for the phase transition. By adjusting T, one can approach the phase transition by tuning 
the mass m of the three-dimensional theory toward zero. The study of second-order (or 
weakly first-order) phase transitions is therefore the study of the infrared behavior of three- 
dimensional theories. 



K 

FIG. 2. The renormalization group flow in Ac/)* theory. 

Now consider this theory in d = 4— e spatial dimensions, instead of 3, where e is finite 
but small compared to 1. The infrared behavior can be studied by following the flow of 
relevant couplings using the renormalization group (RG). For a scalar theory, one has 

/i-^A(/i) = -eA + /3(A). (1.3) 
djj 

The first term, which vanishes in the four-dimensional limit, reflects the trivial classical 
scaling arising from the fact that the interaction / d'^xcf)^ has engineering dimension — e. 
The function /3(A) is the usual four- dimensional /3-function given perturbatively byQ 

/3(A) = ^A2 + 0(A3). (1.4) 

The renormalization group flow is depicted qualitatively in Fig. |^, where the arrows denote 
flows into the infrared {i.e., decreasing /i). There is an infrared- stable fixed point at 

A. = ^'e + 0(e^). (1.5) 

Such a fixed point indicates a second-order phase transition since, if the temperature has 
been adjusted so that m is zero, the behavior of the theory at the fixed point in A will then 
look the same on all distance scales. 



'^Throughout this paper, we use a minimal subtraction renormahzation scheme. As is wefl known, in this 
scheme the /3-function has no additional e-dependence H. 
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Note that the fixed point couphng is small, A <^ 1, provided that e -C 1. For this reason, 
the one-loop result for the /3-function ( |1.4| ) is sufficient to find the fixed-point to leading 
order in e. Renormalization-group improved perturbation theory in A is therefore equivalent 
to an expansion in powers of e. It is for this reason that calculations of the phase transition 
are more tractable in 4— e dimensions, when e is small, than in three dimensions. 

By computing anomalous dimensions of various operators at the fixed point, one can 
extract the critical exponents associated with the phase transition. For instance, one finds 
in a three-loop calculation that the susceptibility exponent 7, defined by x ~ \T — T^y where 
X is the susceptibility, is given by PJIO[] 

7= 1 + 0.167 e + 0.077 e^- 0.049 e3 + 0(e^). (1.6) 

Perturbation expansions are asymptotic in A, and the terms in the expansion start grow- 
ing in magnitude at orders n > 0(1/A). Expansions in e are therefore also asymptotic, with 
terms growing in magnitude at orders n > 0(l/e). What does this imply when one finally 
returns to the three-dimensional theory by sending e ^ 1? If one is lucky, 0(l/e) really 
means something like three or four when e = 1 and the first few terms of the series will 
be useful. If one is unlucky, no terms in the expansion will be useful. Whether or not one 
will be lucky cannot be determined in advance of an actual calculation. From the result 
( |1.6|) , we see that the terms displayed do indeed get smaller for e = 1, and their sum in 
1.195. This compares favorably with results determined by high-order Borel resummation 
techniques |jll[, which give 1.2405 ± 0.0015. (This value is consistent with numerical results 



using lattice methods [0].) So the straightforward use of the e-expansion above appears to 
get within 4% of the true answer — a stunning achievement. 

Another example, which isn't quite as nicely behaved, is the anomalous dimension 77 
of the two-point correlation function, defined by (0(r)0(O)) ~ ^'^-d-v^ Because its leading 
order contribution happens to vanish, the exponent 17 is a small number. One finds P,p!3 



7] = 0.0185 + 0.0187e^ - 0.0083 + 0.0359 + O(e^) . (1.7) 
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If one takes e ^ 1 and adds terms until they start increasing, one obtains rj = 0.0289. The 
correct answer is beheved to be ~ 0.035 [|11],|12[. 



C. Electroweak theory 

We now turn to the apphcation of the e-expansion to electroweak theory.^ We shall ignore 
the Weinberg mixing angle, and concentrate on an SU(2) gauge theory with a single complex 
Higgs doublet. One must first reduce the theory to an effective three-dimensional one, and 
then map out the renormalization group flows in 4— e dimensions. This procedure has been 



analyzed by Ginsparg for general non-Abelian gauge theories ^T^. When constructing the 
three-dimensional theory, one finds that the Aq component of the four-dimensional gauge 
field picks up a Debye screening mass of order gT. The three-dimensional theory then 
consists of the SU(2) Higgs doublet a three-dimensional SU(2) gauge field A{x), and 

a massive scalar Aq{x) that is in the adjoint representation of SU(2). We are interested in 
the infrared behavior of the theory, however, and the massive Aq field will decouple from 
distance scales large compared to 1/gT. One may integrate out its effects as well, leaving a 
three-dimensional SU(2)-Higgs theory of $ and A. 

Now consider instead a 4— e dimensional SU(2)-Higgs theory. The one-loop renormaliza- 
tion group equations have the form 

fi— = -eg^ + Pog'^, (1.8a) 
d/i 

= -eA + (a^^ + S/A + cA^). (1.8b) 

The precise numerical coefficients, as well as the explicit solution, will be given in section 2. 
One finds that the flows have the form shown in Fig. |^. There is no infrared stable fixed 



®A few authors have recently made partial attempts to apply e-expansion techniques to current problems 
in the study of the electroweak phase transition. Ref. examines related problems in the cubic anisotropy 
model, where first-order transitions were originally studied with the e-expansion by Rudnick [p^ . Ref. |l6| 
gives a very rough, heuristic attempt to incorporate renormalization-group results into an order-of-magnitude 
estimate of corrections to the conventional one-loop potential. 
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point, and all trajectories flow to the region A < where the Higgs sector appears to 
become unstable. This instability, and the absence of stable flxed points, suggest that the 
phase transition is always flrst-order in 4— e dimensions [0. 




FIG. 3. The renormalization group flow for an SU(2)-Higgs theory. The dashed hne is the trajectory 
which flows from the initial couphngs {gf, Ai) into the region where \ <^ . 

The last point deserves more careful analysis, which can be provided by recalling the 
previous discussion of the one- loop effective potential. Notice that, although the couplings 
may start at a point [gl^Xi) where the naive loop expansion parameter X/g'^ is not small, 
they eventually flow into the region \ <^ g^. At this point one can apply the familiar analysis 
based on computing the effective potential. As we shall see, the phase transition is indeed 
flrst-order when e is small and its properties can be reliably computed. 

In order to make contact with e = 1, it is useful to consider how Fig. ^ changes with 
e. The one-loop renormalization group equations have the interesting property that the 
trajectories in the {g'^,X) plane are independent of e if g^ and A are rescaled by e, but the 
rate at which those trajectories are traversed is exponentially sensitive to e. Speciflcally, by 
rewriting 




p;2 



A^eA 





the equations become independent of e. 



n'^ = -~X + (ag^ + bg^~X + c~X^). (1.10b) 

Ct/i 

Hence, for these one-loop equations, Fig. ^ does not change with e provided we label the 
axes as g'^/e and A/e. 

We are now in a position to outline our basic approach for using the e-expansion to 
compute quantities associated with the phase transition. This approach is closely related to 
the methods used years ago by Chen, Lubensky, and Nelson to study U(l)-Higgs theory {i.e., 
superconductivity) [|l^ and before that by Rudnick to study the cubic anisotropy model |jl5|.| 



Start by considering the original four dimensional theory and the effective three dimensional 
theory, regulating both with dimensional regularization and renormalizing at the scale fii = 
T. Relate the three dimensional couplings to the four dimensional ones by computing 
perturbatively, in both theories, physical quantities characterized by the momentum scale 
T and equating the results. For example, one can compute the four-point scalar amplitude 
when the external momenta are of order T. One finds a perturbative relation between the 
three- and four-dimensional couplings which is trivial at leading order: 

gl = fiT) + 0{f), (1.11a) 
\, = \(T) + Oif,\''). (1.11b) 

Here, gf and Ai are the three-dimensional couplings at a scale fii = T, and g"^ and A are the 
couplings of the original 4-dimensional theory. The three dimensional couplings gf and Ai 
are the starting points for the renormalization group flow depicted in Fig. ^. Now go from 
three to 4— e dimensions while holding gf/e and Ai/e flxed, so that the relative position in 
Fig. ^ does not change. Finally, use the renormalization group to flow to an equivalent theory 
in the region X <^ g^. The phase transition in this theory can now be studied using a loop 
expansion to compute the effective potential. Note that, for small e, scaling the couplings 
gf and Ai to be 0(e) while holding \i/gf flxed, automatically implies that ^ Ai ^ 1. 



See part II, chapter 4 of Ref. [[19| for a review. 
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In this paper, we shall not concern ourselves with the calculational details of the higher- 
order matching ( |1.11| ) of the three-dimensional couplings with the original four- dimensional 



ones. We shall instead focus on the infrared behavior of the three-dimensional theory as 
predicted by the e-expansion. This is appropriate if the original four-dimensional couplings 
and A are taken to be very small while the ratio X/g"^ is held fixed. So we shall study the 
limit 

^^A^O, with A/^2 fixed. (1.12) 

In practice, this means that we will assume that A ~ ^ 1 but A ^ g^, g^, etc. This 
is appropriate for the real electroweak theory if the Higgs mass is small compared to 1 
TeV but not much smaller than the W-mass. We make this assumption that the couplings 
are small only to simplify the relationship between the four- and three-dimensional theories. 
When working within the three-dimensional theory itself, however, we shall generally present 
formulas that are valid even when gf and Ai are not assumed small. 



D. Outline 

Our presentation in section II begins with the one loop analysis. The explicit solution of 
the one-loop renormalization group flow equations is presented, followed by the evaluation of 
the one-loop effective potential (in 4— e dimensions) when the scalar self-coupling vanishes. 
This information is then used to compute, to one-loop order within the e-expansion, the 
scalar correlation length at the phase transition in both the symmetric and asymmetric 
phases, the baryon violation and bubble nucleation rates at the transition, the latent heat 
of the transition, surface tension, and the free energy difference between symmetric and 
asymmetric phases near the transition. The predictions of the e-expansion for these physical 
quantities are then compared to the results obtained from one-loop perturbation theory 
directly in three dimensions, in the limit where X <^ g^. 

Section III contains two loop analysis. The two-loop renormalization group equations are 
presented, along with their explicit solutions, followed by the two-loop effective potential for 
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A = 0. Two-loop corrections to the scalar correlation length and latent heat are computed. 
We again compare results with direct three-dimensional perturbation theory in the \ <^ 
limit. More significantly, by examining the relative size of the two- loop corrections, we can 
test whether the e-expansion is well-behaved in the case A > (7^ of moderate-to-heavy Higgs 
mass, where conventional perturbation theory fails. 

Section IV examines the e-expansion when one increases the number of scalar fields, A^. 
Near four dimensions, when the number of charged scalar fields is sufficiently large {N > Nc), 
two additional renormalization group fixed points appear: an infrared stable fixed point and 
a tricritical point. We compute the first two terms in the e-expansion of the critical value 
A^c at which these fixed points first appear. We also examine the N ^ 00 limit of the 
slope X/g"^ of the tricritical line separating theories with first and second order transitions, 
and compare the e-expansion result with that obtained directly in three-dimensions using 
large- methods. 

The overall interpretation of our results is summarized in the conclusion. Various details 
of two-loop calculations are relegated to appendices. 



E. Possible Misgivings 

Before embarking on explicit calculations using the e-expansion, it is worth considering 
whether the endeavor is doomed before it begins. In particular, the e-expansion predicts 
that the phase transition is always first order, even when Xi/gf is arbitrarily large. Some 
authors have dismissed this result based on the apparent failure of analogous arguments for 
superconductivity.^ At large distances, simple superconductors are described by Landau- 
Ginzburg theory, which is just U(l) gauge theory with a single, charged, complex scalar field. 
The renormalization group equations in 4— e dimensions have the same form as Eq. ( |1.8|) 
with jSo positive. The fiow is shown in Fig. HI. As observed by Halperin et ai, pT|, the e- 
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For a nice review of such arguments, see Ref. EC 
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expansion predicts the superconducting phase transition to be always first order, even deep 
in the Type II regime, Ai ^ gf. 
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FIG. 4. The rcnormalization group flow for a U(l)-Higgs theory with N < Nc 



The predicted strength of the transition is too small to measure in superconductors, 
but detailed studies have been made of the critical properties of smectic-A liquid crystals, 
which are argued to be in the same universality class as superconductors |Q . The predicted 
first-order transition is often not observed. Motivated by this discrepancy, Dasgupta and 
Halperin made a lattice study of the superconducting transition. They used a U(l) 
lattice gauge theory in three dimensions with a fixed-length Higgs field; that is, their bare 
theory had Xi (yo with m^/ Ai fixed. Simulating a bare theory that had gf = 5, they found 
no numerical evidence of a first-order phase transition. 

Their numerical results are not ironclad. The height of the specific heat peak was mea- 
sured on lattices of size 3'^, 5^, 10^, and 15'^. In a first-order phase transition, the height 
should grow like the volume V in the large- volume limit; for a second-order transition, it 
should grow like V^"''^ (or InV if the specific heat exponent a ~ 0) p3[. Results on these 
four lattice sizes did not support a linear relationship between the height and V. However 
it is not clear if these simulations reached the large volume limit. As a rough guide, one 
may consider the predictions of the e-expansion. If we use the one-loop flows and set e = 1, 
we find that the scale change necessary for the initial theory with Xi = oo and gf = 5 
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to run to the region A < is roughly 9. (Detailed formulas are given in section 2.) This 
suggests that the width of domain walls, or the minimal size of critical bubbles nucleated in 
a first-order phase transition, will be of order 9 lattice units. Therefore, the linear growth 
of the specific- heat peak with V need not begin until the lattice size is larger than 0(9^). 
Given that the lattices of |^ range in size from 3^ to 15^, the reliability of the conclusions 
is clearly sensitive to whether the condition V ^ 0(9'^) really means something like V > 4^ 
versus something like V > 18^. 

Nevertheless, given the other evidence from liquid crystals, let us put this concern aside 
and hypothesize that the phase transition studied in pS) is in fact second order. The 
presence of a continuous transition in a fixed-length U(l)-Higgs theory would imply that 
the e-expansion does a terrible job describing the flow of theories starting near Ai = oo. 
These theories must flow to some infrared- stable fixed point not seen within the e-expansion. 
However, as noted after Eq. ( |1.12|) , in the application to electroweak theory one is interested 
not in theories that start at strong coupling with Ai = oo, but in perturbative theories 
that start near the unstable Gaussian fixed point at X = = 0. The latter theories need 
not lie within the domain of attraction of the putative fixed point, and might instead be 
correctly described by the e-expansion. Within the e-expansion, all theories which start 
near X = = stay inside the shaded region of Fig. ^ Recall that at leading order the 
e-expansion is based on renormalization-group improved one-loop perturbation theory. One 
way to view the success of the e-expansion in pure A0^ theory is to note that, even though 
the fixed point coupling A is formally 0(1) when e — 1, the coupling is not such a large 
number that higher-order corrections in A cause dramatic changes. If we are similarly lucky, 
the couplings in the shaded region of Fig. ^ will not be so large when e — 1 as to render the 
e-expansion useless. 
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FIG. 5. The renormalization group flow for a U(l)-Higgs theory with > A'c- Point A is an in- 
frared-stable fixed point, while i? is a tricritical point. The dashed line is the tricritical line which separates 
theories with first and second order transitions. 

An important general point is that phase transitions are more hkely to be second order in 
lower dimensions due to the increased importance of long-distance fluctuations. A prediction 
of a second-order phase transition in 4— e dimensions may be more robust as e — > 1 than 
a prediction of a first-order phase transition.^ A useful setting for discussing this point 
is provided by the generalization of the single scalar U(l)-Higgs theory to a theory of 
charged scalars with a global U(A^) symmetry. For A^ very large, one finds that the flow 
looks like Fig. |^ instead of Fig. ^. An infrared-stable fixed point and a tricritical point 
have appeared. Theories which start to the right of the dashed line have second-order 
phase transitions; those to the left undergo first-order transitions. Near four dimensions, 
the transition between Fig. ^ and Fig. ^ occurs at A^c(4) = 183.0. As the dimension d is 
decreased and second-order transitions become more likely, the critical number of fields, 
Nc{d), should decrease. We would like to know whether A"c(3) < 1. 

In section IV, we compute the first sub-leading correction to A"c within the e-expansion. 
This is only a computation of the slope of Nc{d) near d = 4. However, for both U(l) 



"'^^For example, as an alternative to the d = 4— e expansion, some work has been done using ad — 2+e 
expansion |2j]. These studies find second-order phase transitions, although they are restricted to considering 
theories with Ai = oo. 
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and SU(2), when extrapolated to e = 1, we find that the 0(e) term is larger than the 
leading order term and of opposite sign, so that A^c naively extrapolates to A^c(3) < 1. But 
one cannot trust this result because the series is clearly badly behaved. This in itself is 
somewhat discouraging because A^^c is revealed to be an example of a physical quantity for 
which the e-expansion is not very well behaved. 

We do not believe that there is any compelling evidence one way or the other as to 
whether the SU(2) phase transition, with one Higgs doublet, is first or second order for large 
Higgs mass. This is an issue to be resolved by lattice simulations or other techniques. If it is 
second order for large or moderate Xi/gf, then the e-expansion is doomed. Our philosophy 
will be to assume that the transition is indeed first-order and see what we can extract from 
the e-expansion. 

Before closing this section, we note one other hurdle for the e-expansion posed by recent 
literature. One can solve the U(l) and SU(2) theories directly in three dimensions in the 
limit ^ oo. Alternatively, one can compute to some finite order in the e-expansion, set 
e = 1, and only then take the N oo limit of the result. A comparison of the two results 
gives a test of the e-expansion. In particular, one can compute the value of X/g"^ for the line 
that divides first- from second-order transitions. This has been investigated in the large- A^ 
expansion by Jain et al. |25|. Their result is completely different from that suggested by the 
e-expansion, differing even in the power of A^. In section IV, we discuss the flaw in their 
analysis and redo their calculation. We find that the critical value of \/ g'^ scales like 1/A^ 
as A^ ^ oo and that this behavior is in agreement with the e-expansion. 



F. The magnetic mass and non-perturbative physics 

It is well known that unbroken non-Abelian gauge theories at high temperature (such as 
SU(2) electroweak theory in the symmetric = phase) suffer a breakdown of perturbation 
theory beyond the first few orders. The loop expansion parameter is not just large; it is 
infrared divergent. This reflects the fundamentally non-perturbative nature of long distance 
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physics in these theories. In this section, we discuss how this affects the e-expansion. 




FIG. 6. An example of an infrared divergent contribution to the transverse gauge field mass in three 
dimensions. 

An easy way to see the problem is to study the effective three-dimensional theory and 
attempt to compute the transverse gauge field mass in perturbation theory. Gauge invariance 
and rotational invariance imply that the self-energy of the gauge field A should have the 
form 

n.,(p) = (p%-p.p,)/(p'). (1.13) 

Unless f{p'^) singular as p'^ 0, there will be no mass at any order in perturbation theory. 
However, fig. ^ shows a self energy graph which is infrared divergent. This divergence can 
be seen using simple three-dimensional power counting which shows this diagram gives the 
logarithmically-divergent contribution 

at zero external momentum. Higher order diagrams have increasingly severe divergences. 
Non-perturbative fluctuations produce a finite correlation length which serves to cut-off 
these divergences. The inverse correlation length (or magnetic mass) is of order gfT (or 
(7^yu). If one now computes a three-loop diagram using this scale as an infrared cut-off, one 
finds a contribution of order 



d^ 

glT 

This is the same size as the two-loop diagram, as are all higher-loop contributions! Pertur- 
bation theory is therefore inadequate to treat infrared physics at scales of order g\T. 
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What happens to this problem in 4— e dimensions? The two-loop diagram of fig. |^ is no 
longer infrared divergent when e < 1, and infrared problems do not arise until higher order. 
In particular, an n-loop contribution has the form 

(^vr/^, (1.16) 

so that infrared divergences first appear at order n = 2/e. Perturbation theory is useless 
beyond this order, which means that no quantity can be computed perturbatively with 
relative error smaller than gf" = g"^^^. Since is 0(e) in our procedure for the e-expansion, 
this relative error is [0(e)]^/^. But this is just the same order as the uncertainty intrinsic 
to the e-expansion in any theory due to the fact that the expansion in e is an asymptotic 
series; terms of order e^ in the series start getting larger beyond order ~ 1/e. So the non- 
perturbative physics related to the mass gap in three-dimensional non-Abelian gauge theories 
may not pose any problem for the e-expansion not already present from the asymptotic 
nature of the expansion. 

Some physical quantities we will examine are infrared-safe as they depend only on physics 
in the asymmetric phase (where all particles have perturbative masses), or do not probe the 
symmetric phase over long distances.^ Other quantities, such as the scalar correlation 
length in the symmetric phase are, in three dimensions, infrared safe at lowest order but not 
at higher orders. How well the e-expansion will work, after extrapolating e 1, for such 
quantities remains unclear. This issue will be discussed further in section III C. We will, 
however, largely focus on infrared-safe quantities when computing beyond leading order. 

II. ONE-LOOP ANALYSIS 



The sphaleron mass (or baryon violation rate) in the asymmetric phase is an example of the latter case. 
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A. Renormalization group flow 



We are now ready to look in detail at the solutions to the one-loop renormalization- 
group equations. We shall quote results for both U(l) gauge theory with N charged U(A^)- 
symmetric scalars, and SU(2) gauge theory with N/2 U(A^/2) -symmetric scalar doublets.^ 
will always denote the number of complex scalar degrees of freedom. The Lagrangians 
are[^ 

^U(l) _ 1 -~ - ■ - ^ 2 



-E - 2 

3U{2) _ 1 
-E - 2 



(9-^/i^/V)$ +>'|<fp + iT/i^A|$r + iF2^, (2.1a) 
£SU{2)^i (9-^^^/VV,/2)$|' + lm2|$|2 + l;x^A|$r + i (f;,)'. (2.1b) 



We will express potentials in terms of the scalar field modulus = |$| = {Y.a^aY^'^i 
terms of which the classical potential takes the canonical form, 

= > V + i /^^A 0^ (2.2) 

We choose couplings A and g"^ to be dimensionless in any spacetime dimension and therefore 
need to include explicit powers of the renormalization point fi in the interaction terms. 

The renormalization group will be used to run from a renormalization point /ii = T to a 
final scale /i. Our notation will later be more compact if we write /i = /Ui/s, where s starts 
at 1 and grows to oo as one flows into the infrared. In place of the gauge couplings in ( p.l| ) 
with standard normalizations, it will also be convenient to use the "charge" , 

, 9, U(l) 

q ^ (2.3) 
9/2, SU(2) . 



-'^^More specifically, in tlic SU(2) theory we require Sp(A^)r custodial symmetry, which is the U(iV/2) flavor- 
symmetric generalization of the usual SU(2)r ~ Sp(2)r custodial symmetry of the one-doublet case. The 
complete global symmetry SU(2)l x Sp(A^)r implies that the scalar potential can only be a function of a 
single variable, J2a \^a\'^j so that the potential is actually an 0(2A^) invariant function of <!>. 



^^We use an unconventional normalization in the Lagrangians (2.1); a more conventional choice is obtained 
by rescaling the field $ by V2. 
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The mass of the gauge boson is then q/ji^cj) in both theories. The introduction of the charge q 
will allow us to write common expressions for most one-loop results in either gauge theory. 
The one-loop RG equations take the form 

s^ = eq'-/3oq\ (2.4a) 
s^^eX-{aq^ + bq^X + cX^), (2.4b) 

where the constants are given by 

^° " 3(47r)2 ' ^"(4^0^' (4^0^' 3(47r)2 ' ^^"^^ 

where 

[1, U(l); 

% ^ (2.6) 
1 3, SU(2), 

is the number of gauge bosons. C2 is the quadratic Casimir which, for our two gauge groups 
of interest, equals ng — 1. We shall later also need the renormalization group equation for 
the mass m^. 



ds 



-{loq' + hX)m\ (2.7) 



where 



6ng , _ 2(Ar + l) 

^° - " (4^0^ ' - ■ ^^-^^ 

We shall not need the renormalization group equation for the field (f). 
The solution of the RG equation for q"^ is 



q\s)^s^qi/ [l + ^{s^-l)qi\. (2.9) 

The solution for A is most easily found by changing variables to the ratio x = X/q^ whose 
RG equation, 

= + (5 _ p^)^ + cx% (2.10) 

ds 
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has the solutionis 



^is) = — {(3o-b+VA tane{s)) 



where 



1 + —{s -l)q^ 



_i ^ 20X1 - {l3o - by 



and 



a = Tan 



A = 4ac - {Po - hf . 



(2.11) 

(2.12) 
(2.13) 

(2.14) 



This is the form of the solution most convenient for the cases of Fig. ^ or Fig. ^, which occur 
when A > and no stable fixed point exists. When A < 0, two new fixed points appear as 
shown in Fig. ^, and Eq. ( p.ll ) may be put in the equivalent form: 

wis) — 1 



X[S) 



1 



wis) + 1 



where 



and 



w{s) = 



Xi — x^ 
X+-X, 



1 + — (s -l)gi 



-A/A) 



(2.15) 



(2.16) 



(2.17) 



The lines x = x i and respectively through the stable fixed point and the 

tricritical point shown in Fig. ^. 

Return now to the solution ( |2.11D describing the purely first-order case with which we 
are mostly concerned. We want to flow to the region A <^ g^, where we will use perturbation 
theory. For the sake of simplifying formulas and derivations, it is particularly convenient to 
run precisely to A = 0. When A(s) = 0, the values of the scale factor s and the charge 
extracted from Eq. (|2.11|) areQ 



^^This solution was previously derived for the U(l) case in jist - 
In Ref. it was incorrectly asserted that the strength of the phase transition is exponentially small 
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and 



l + 4fexp 



where we have defined 



Wo 



2Po 



2 \ x^^/A ^ 



2 V x^VA , 



(2.18) 



(2.19) 



^0^ 



/So' 



(2.20) 



The Xi — > oo hmits of these formulae were used to derive the estimate of the scahng factor 
appearing in our discussion of the lattice simulations of Dasgupta and Halperin in section 
^E| . By taking both the xi ^ oo and gi ^ limits, one can extract the maximum value of 
q in the shaded region of Fig. H, 



% 1 - exp 



^ + Tan-(^ 
2 V a/A . 



(2.21) 



B. The one- loop potential 

Having run to A ^ g^, and in particular A = 0, we now want to use ordinary perturbation 
theory to compute various quantities describing the phase transition. The tree-level Higgs 
potential is 



yW(0) 



1 + -M + hX- 27^) 
e 



A + -(ag^ + bq'^X + cA^ - 4A7<i) 
e 



(2.22) 



when qi 1 and Ai/qJ > 1. The assertion was based on examining the form of ( 2.10|) and noting that for 
X ^ 1 the right-hand side is of order q^. To run a; to requires a change in x of order 1 and hence, roughly, 
a change in Ins of order 1/q^. Ref. fl^ then asserted that his must be order l/qf, implying that s must run 
exponentially far before x becomes small. This is incorrect because q{s) is not a constant. From Eq. ( ^.9D , 
q{s) runs to be order 1 when s*^ is order l/qf, and at that point x can quickly run to zero. But s = 0{qi ^^'^) 
is not exponentially large when e ^ 1. 
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where we have included the one-loop counter-terms for the Minimal Subtraction (MS) scheme 
and 70 = 70(A,g^) is the one-loop anomalous dimension of 0. Working in Landau gauge, 
the one-loop potential is 

V^(i)(0) = V^(O)(0) +X(rn2+iA02) + {2N-l)X{m'+\\4>'') + {'i-t)X{q^^) , (2.23) 
where it is notationally convenient to introduce the dimension one field 

= /i^/^- (2.24) 

The basic one-loop contribution I is given by the (dimensionally regularized) integral 

^ r d'^k ^ , r(-2 + e/2) , , , 

The arguments of the three X functions appearing in Eq. ( p.23|) are the masses of, re- 



spectively, the Higgs boson, the unphysical Goldstone boson, and the vector boson in the 
background of (p. 

There are a few technical points to mention before proceeding. First, we have ignored 
any 0-independent constants in the effective potential. Such constants are not relevant to 
any of the quantities we shall be considering; including the constant term would require 
that we study its renormalization-group flow as well. Second, we emphasize again that 
our convention is d = 4— e and not d = A — 2e. Finally, we have used MS regularization 
rather than, say, MS renormalization. Including finite counter-terms to implement MS 
regularization would not, of course, change the relationship between physical quantities. 
For example, an expansion of critical exponents of a second-order phase transition in terms 
of e = 4 — c?, which is "physical," would be independent of whether intermediate calculations 
were done in the MS or MS scheme. However, in our application final results will also 
depend on Ai and qf, whose definition are scheme dependent. But the difference between 
the MS and MS definitions of Ai and is order 0(A^, g^), and hence can be ignored under 
our working assumption ( |1.12| ) that Ai and qj are small (with Ai ^ qf). 



Now substitute A = into the above potential. Ignoring 0-independent constants, only 
the vector contribution to the one-loop potential survives and one finds 
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1 - Me) 



where 



Me)^e(l--)r(-2 + ie)e^W2^ 



(2.26) 



(2.27) 



and where it's convenient to introduce the MS scale 



47r 



Some useful limits of h{e) are 



3 tt' 



(2.28) 



h{e) = 1 + — e+ — + 



12 \ 16 4J 



Oi 



and 



Ml) 



167r /eTE 



9 V 47r 



(2.29) 



(2.30) 



Expanded to leading order in e, the one-loop potential is 



In 



fx 



5 

12 



0(6). 



(2.31) 



However, for the sake of later tests of the e-expansion, it will be convenient to keep the 
potential in the more general form of Eq. ( p.26|) . Note that substituting e = 1 into (p.26|) 
produces the cubic 0^ term familiar from earlier studies of the electroweak phase transition 
employing the one-loop effective potential |jl]J|j2^. 

Recall that the three-dimensional mass m is adjusted by changing the temperature in the 
original, four- dimensional theory. At the phase transition, the potential has two degenerate 
minima, so that V{(j)c) = ^(0) and {d/d(j))V{(j)c) = where 0c is the non-zero Higgs con- 
densate at the asymmetric minimum. Applying these conditions to the one-loop potential 
(g) yields 



2 a^V. 



and 



= -/(e) 



(2.32) 



(2.33) 
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where 

1/6 



e-i/12 



(2.34) 



It should be emphasized that these are the values of m'^{fi) and </'(yu) at the scale fi where 
A(yu) = 0, and not at the original scale fii = T. 

Before using these results to examine various physical quantities, we should make a 
few comments about our choice of scale so that A(s) = 0. First, recall that the vector 
contributions to the potential are the important ones, and so the most important physical 
scale which affects the phase transition is the vector mass M ~ q(j). One should therefore 
choose /i ~ M in order to avoid the appearance of what in four dimensions would be large 
logarithms ln(M/ /i) in the perturbation expansion and in 4— e dimensions are large factors of 
[(M//i)^ — l]/e. But, a posteriori, we see from ( p.33| ) that this is precisely what we have done. 
This is an alternative phrasing for our previous criteria that perturbation theory is good when 
A/g^ is small. Slightly different choices of /i ~ M (for example, choosing A(s) = 0.1 q^{s)) 
will give the same results for physical quantities, order by order in perturbation theory. 

But what would happen if we ran just a tiny bit further until A(s) < while keeping 
A/g^ small? The classical potential ( |2.2| ) is then unbounded below! This disaster is illusory, 
however, because both the bare potential ( p.22|) with one-loop counter-terms and the one- 



loop potential (|2.23| ) remain bounded below. It must be remembered that A in the MS 
scheme is not itself a physical quantity, and there is no reason why it can't be negative for 
a sensible theory. This is why the running of A to negative values was only suggestive of a 
first-order transition for small e. Only by computing the effective potential, when |A| <^ g^, 
can one be sure. 



C. Applications 

We are now in a position to put everything together and compute various quantities of 
physical interest. The most interesting qualitative conclusions will follow from our first two 
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examples, the scalar correlation length and the rate of baryon number violation. We will also 
summarize results for a number of other quantities; these results will be useful when testing 
the e-expansion against unimproved perturbation theory {i.e., conventional perturbation 
theory without the renormalization group) at the end of this section. For the sake of these 
later tests, it will be helpful to have on hand full one-loop results in any dimension, rather 
than just the first term or two in the expansion around e = 0. 



1. Scalar correlation length, ^ 

As a first example, consider the scalar correlation length ^sym = l/'^c in the symmetric 
= phase at the critical temperature. Using the one-loop result ( |2.32D for the critical 
mass, and the definitions /i = /ii/s = T/s, we have 



= I (l - i) {1 + Oiq\s))}, (2.35) 

where q{s) and s are given by Eqs. ( p.l8| ) and ( |2.19| ). Similarly, the correlation length at the 
critical temperature in the asymmetric phase where = (/)c is obtained from the curvature 
of the one- loop potential at (j)c- 



Uym = ^ {1 + OiqXs))}. (2.36) 

Recall that our implementation of the e-expansion involves holding qf/e and Ai/e fixed. 
Hence qf, q'^{s), q^ and Ai are all of order e. From Eq. ( p^.l9D , it follows that s"^ is order e° 



and so the correlation length in either phase depends on e as 

^^e-i/2eO(iA). (2.37) 

For small e, the dominant effect is clearly the exponential dependence of s on e. The e- 
dependence is not a simple power series expansion like we reviewed for critical exponents in 
section It can be made similar by taking the logarithm: 

+ 0(e) (2.38) 



Inf s'') 

In(easymT) = — ^ - ^ lu 6 + 



^ \87rag^(s) / 



0(e-i)+0(lne) + 0(e°) + 
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Note that terms of 0(e°) in this formula which come from the computation of the one- loop 
potential are of the same order as the s«6-leading corrections to s'^ in the first term. So, 
by the philosophy of the e-expansion, we should not keep the 0{e^) terms unless we also 
compute the corrections to s"^ using the two-loop renormalization group. We shall do so in 
section 3. Note that the leading 0(l/e) term in ( |2.3(j| ) gives nothing more than ^ ~ s/T. 

Because this lowest-order calculation does not determine the overall normalization of the 
correlation length {e.g., it cannot distinguish between ^ ~ s/qoT and ^ ~ s/(47r^go^)), it 
is not terribly useful for practical applications. However, one interesting comparison that 
is insensitive to an overall (e-independent) multiplicative factor is the ratio of the above 
result to what would have been found with straight one-loop perturbation theory in 4— e 
dimensions without using the renormalization group. The exponential dependence on 1/e 
of the RG unimproved result can be extracted by examining the one-loop potential (p.23|) 
or, equivalently, by noting that the perturbative solution to the RG equation ( p.lO|) is 



x{s) = xi - qf{a + [b - (3o)xi + cxf) ( ) + O 



(2.39) 



If only the first two terms in this expansion are retained, then the solution for x{s) = gives 

= qf{a+{b-f3o)x, + cxl)- ^^'^^^ 




•iS. 0.6 

OJ 

to 

0.4 
0.2 

^ ' 
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FIG. 7. The ratio of the RG-improved to unimproved values of the scale factor s'^ at one loop order for 
electroweak theory. 
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In Fig. ^ we plot the ratio of s'^/s^^j.^ versus Xi/qf for the SU(2) theory with a sin- 
gle Higgs doublet (A^ = 2). is the RG-improved result of Eq. ( |2.19| ). We have fixed 
qf/e to be the value (0.32)^ of electroweak theory and translated Xi/qf into the zero- 
temperature Higgs mass mh(0) using the tree-level relationship m1{0) = (AXi/Sgf) M^(0) 
with Mw(0) = 80 GeV. Fig. |^ suggests that, for light and moderate-mass Higgs bosons (up 
to roughly 130 GeV), the phase transition will be stronger than implied by the unimproved 
one-loop calculation. For heavier Higgs bosons, the transition will be weaker. The increased 
strength for light mass Higgs is qualitatively consistent with the results of explicit two-loop 
calculations in Refs. 

Typically, most quantities with positive mass dimension, such as the latent heat, the 
Higgs expectation value (in a particular gauge), or the height of the free energy density 
barrier between the two vacua, will scale with a positive power of s. At leading order in the 
e-expansion, when compared to unimproved one-loop calculations, they will therefore have 
a behavior similar to that shown in Fig. |^ We shall examine these quantities in more detail 
later. First, however, we will examine the potentially more subtle behavior of the rate of 
baryon number violation. 

2. Baryon violation rate, Fb 

The original motivation for this study was to understand whether baryon number viola- 
tion is sufficiently slow in the asymmetric phase to be compatible with electroweak baryo- 
genesis. For simplicity, we shall study this rate exactly at the critical temperature where 
the two vacua are degenerate. In the case of the real electroweak phase transition, there is 
actually a small amount of supercooling below this temperature before the phase transition 
completes. 

Since the rate of baryon number non-conservation per unit volume has positive mass 
dimension, it will naively scale with a positive power of s just like the latent heat or the 
Higgs expectation. By itself, this suggests that when the phase transition is stronger than 



30 



predicted by an unimproved one-loop calculation, then the baryon violation rate will also 
be faster. However, this argument is overly simplistic. The scale factor s is exponentially 
sensitive to 1/e, which is why the dependence on s dominates the small e dependence for 
most quantities. As we shall discuss, the baryon violation rate has additional exponential 
sensitivity to Since is 0(e), this provides an additional exponential dependence 

of the rate on 1/e that must be properly accounted for when analyzing the rate to leading 
order in the e-expansion. 

The rate of baryon number violation, Fb, is determined by the action of the electroweak 
sphaleron solution . To compute this rate within the e-expansion, we need the action 

of the sphaleron solution in the 4— e dimensional theory. This poses a problem because the 
sphaleron is intrinsically a three-dimensional object, intimately related to topologically non- 
trivial mappings of SU(2) S^. There is no obvious generalization of the three-dimensional 
solution to a rotationally-symmetric solution in 4— e dimensions. On a concrete level, the 
problem can be highlighted by considering the standard form of the sphaleron solution in 
three dimensions 



A{r) = -F{r)r X a, <^{r) = H{r) r ■ a (^^^^ . (2.41) 

There is no natural generalization of f ■ a and f x a away from three dimensions. 

In dimensions larger than three, we may instead consider solutions which have non- 
trivial structure in three dimensions and are translationally invariant in the remaining d—3 
dimensions:]^ 

A{x) = A{xi,X2, X3), $(x) = $(xi, X2, X3). (2.42) 

Instead of computing the sphaleron action Se, we shall compute the sphaleron action per 
unit ((i— 3)-volume, Se/C'^^^, where we imagine putting the system in a very large box of 
size £. 



This is reminiscent of the contortions one must go through to define 75 in dimensional regularization, 
where four dimensions must be treated differently from the others. 
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The action of our sphaleron is determined by the appropriate stationary point of the 
functional 

5e[$, A] = C'-^ jd'x[\ m? + i Fl + . (2.43) 

As before, we shall first run to A(s) =0 so that perturbative corrections will be small when 
e is small. Note that at this point the one-loop potential, as given by Eqs. ( p.26|) and ( p. 321) , 



depends on the coupling q{s) only in the combination qip. We can make explicit all the 
parameter dependence at Tc by writing V^^^cf)) = iJ,'^~^v{q(j)/ ij}""^'^). Now rescale all fields 
and dimensions, 

A^^ A, $, X^-, (2.44) 

q q ^^ 

to produce 

Se = ^^^J ' / rf'x |i |(V, - tT-A,)^\^ + i(Vii, - VjA, + Ax A,y + g2^;(0)| . 

(2.45) 

Recall that q^{s) is 0(e). From this form of the action it is apparent that the details of the 
one-loop potential f (0) do not affect the sphaleron action at leading order in e. However, 
the solution is already known in the case where the strength of the potential is negligible: 
it is the standard sphaleron solution for the classical potential A(0^— t;^)^ when A 0. 
The action of the three dimensional sphaleron solution in this limit is Se = 3.0405 Mw/q!w 
28| , |29[| . In our case, this translates to 



Se = 3.0405 {fiCy-' — (1 + 0(e)) 



3.0405 ifiCy-' (1 + 0(e)) 

23.344 



g2 



(/x/:)^-(l + 0(e)). (2.46) 



The only importance of the potential is to determine the asymptotic Higgs expectation value 
0c through Eq. (|2:33| ). 
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In three spatial dimensions, the rate of baryon number violation per unit correlation time, 
^Fb, is roughly (V/.^'^) exp(— 5e) where V is the volume of space and C, is the correlation 
length. (The distinction between the scalar and vector correlation lengths, and the presence 
of additional prefactors to the exponential, will not matter at leading order in e.) For the 
moment, let us blithely ignore the {^CY~'^ factor in Eq. ( |2.46| ) and recall that ^ ~ s/T. 
Then the leading-order result in the e-expansion is 



V 



T exp 



23.344 



■In(s^) + 0(e° 



(2.47) 



with q^{s) and given by Eqs. (|2.18| ) and ( p.l9|) . Consider comparing this result to a 
straight perturbative calculation that was not improved by the renormalization group. We 
have already discussed this comparison for the explicit dependence on s. For g^(s), the 
corresponding unimproved result is Spg^"?!- A plot of {s) / {s''^^^^^q\) versus \i/q\ is given in 
Fig. p. Inclusion of the renormalization group flow always increases which drives the 

rate Fb larger. In Fig. ^, we put everything together and plot the ratio of the exponent in 
Eq. ( p.47| ) to the result from unimproved perturbation theory. 




50 75 100 125 150 175 200 
mh(0) [GeV] 

FIG. 8. The ratio of the RG-improved to the unimproved value of the gauge couphng q{s) at one loop 
order for electroweak theory. 
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FIG. 9. The ratio of the RG- improved to the unimproved value of the logarithm of the baryon violation 
rate at one loop order. 

Our leading-order result of Fig. ^ suggests that the baryon violation rate is always larger 
than that deduced from one-loop perturbation theory. If this is a correct qualitative descrip- 
tion of what happens when e — 1, then the one-loop bound of mh(0) < 35 GeV derived by 
Dine et al. P|, is indeed an upper bound on the Higgs mass for electroweak baryogenesis in 
the minimal standard model. 

Having reached this punch line, we must return to an annoying technical point concerning 
the factor in Eq. (|2.46|) which we ignored. Remember that fi = T/s, so this factor 

is exponentially sensitive to 1/e but equals unity at e = 1. If we retain the factor, it changes 
the relative importance of the iSe term and the Ins term in the exponent of ( p.47|) in the 



e-expansion. Also, for e < 1, these factors have different dependencies on the size C of 
the system in the extra 1— e dimensions, and so are not comparable for C — > oo. When one 
works to higher orders in the e-expansion some corrections to In Fb/V will be proportional to 
{lj,£)^~'^, and serve to modify the sphaleron action, while other corrections are independent 
of C and modify the zero-mode prefactor. This suggests a natural approach for defining 
the e-expansion of the full baryon violation rate: separately extrapolate to e = 1 the terms 
independent of C, and the (coefficient of) terms proportional to Combine the two 

contributions to IuFb only at e = 1 where {fiCy~^ may be replaced by unity. This procedure. 
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to lowest order, yields Eq. (|2.47|) (evaluated at e = 1) and the results shown in Fig. |.[^ 

Performing a complete next-to-leading order calculation remains a problem for the fu- 
ture. However, as an example of corrections to the sphaleron action, consider the one-loop 
correction due to quantum fluctuations in modes other than the translational zero-modes. 
The prefactor to exp(— iSe) depends on the determinant of the curvature operator fi^pj^ de- 
scribing small fluctuations about the sphaleron solution. The non-zero mode contribution 
to the prefactor has the formj^ 



^ /DetU^V'^' 



iu -\- Avac 



f2.48) 



where ^"^^^ corresponding curvature operator around the ground state, the determi- 

nants are taken over fluctuations in 4— e dimensions, u is the momentum of a fluctuation in 
the 1— e transverse dimensions in which the sphaleron is static, and {A*-*-*} are the eigenvalues 
of Q"^ in three dimensions. The prime on the determinant indicates the omission of the zero 
eigenvalues, which can't be treated by Gaussian integration. Since the only scale in the 
problem is /x, this prefactor must have the form of 

P~exp{(^/:)i-7(g,e)}. (2.49) 



The exponent has the same (/x£)^~^ factor as does the sphaleron action in ( ^.46] ), and so 
InP generates an O(g^) = 0(e) correction to the classical action. 

Another possible approach for defining the e-expansion of the baryon violation rate is to 
choose the length C of the extra 1— e dimensions to be order This allows the zero-mode 
prefactors to be treated on the same footing as the action and other prefactors, but at the 
price of introducing an unphysical dimensionless parameter, Any calculation to finite 



^^The qualitative behavior of the ratio plotted in Fig. ^ is, in fact, insensitive to the presence or absence 
of the 31ns zero-mode contribution. 



-'^^For a discussion of these determinants in the three dimensional case, see, for example, Refs. |Q and 31 
and references therein. 

^°This approach also has an advantage of allowing one to generalize the sphaleron determinant to 4— e 
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order in e would depend on the precise choice of in direct analogy to the dependence 
of calculations on the renormalization scale in conventional perturbation theory. As in 
conventional perturbation theory, the result should become less and less sensitive to the 
precise choice of /i£ when one computes to higher and higher orders in e. Which approach 
will yield the best systematic approximation for the baryon violation rate is not yet clear. 
We emphasize that the complications discussed here are specific to the baryon violation rate 
and will not be present for other quantities discussed below. 



3. Latent heat, AQ 

The latent heat (per unit volume) of the phase transition is defined as the difference in 
energy density of the two coexisting phases at T^., 

JAF 



AQ = AE 



= TAS 

T=T, 



-T- 



(2.50) 

T=rc 



Here, AE, AS, and AF denote respectively the energy, entropy, or free energy densities in 
the symmetric (0 = 0) phase minus that in the asymmetric (0 = (pc) phase. The free energy 
density F is simply TV^fj, where V^s is the effective potential in the three-dimensional theory. 
Understanding the dominant sources of temperature dependence in Ves requires returning 
to the relationship between the effective three-dimensional theory and the original four- 
dimensional theory at high temperature. There is explicit dependence from the initial choice 
of renormalization scale fii in the three-dimensional theory as equal to the temperature T. 
There is also implicit dependence from the variation of the effective three-dimensional mass 
mi = m(/ii) with T, as was shown explicitly for pure scalar theory in ( |1.2| ). In gauge 
theories, this relationship has the form 



m^ 



+ (c, + CA X)T^ + Oiq'T', \"^T^), (2.51) 



dimensions while retaining only a single negative mode. In contrast, taking C ^ oo generates a continuum 
of negative eigenmodes which contribute an imaginary part to InT' proportional to i(/i£)^^'^7r/2. This 
produces the expected prefactor of i, reflecting a single unstable mode, only when e ^ 1. 
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where Cg and cx are constants. Other sources of temperature dependence, such as the 
running of the original four-dimensional couplings or the coefficients of higher-dimensional 
interactions in the three-dimensional theory, will not contribute at leading order. Now, on 
dimensional grounds, the three-dimensional result for AVes must have the form fif f{ml/ fif). 
Combining this with Eq. ( [^.51D , plus the condition that AV^gf(Tc) = (by definition of Tc), 
allows one to rewrite the latent heat as 



AQ 



-T 



2(0 dm\ d 



cfT 



fj.^=T=Tc 



fJ'l- AVeff 



2 

mf 



(2.52) 



neglecting corrections of order gi or V^i. The mass mi should be held fixed in the above 
derivative. 

We can now turn to the computation of the (rescaled) latent heat AQ/u'^T in 4— e 
dimensions. We have previously determined m(/i) at the phase transition (Eq. ( p.32|) ). The 
mass mi is related to m(/i) by the renormalization group equation ( p.7|) , whose solution is 



m^(/i) = ml P(s^ 



where at one loop 



Pis' 



exp 



1 
e Ji 



iiKs'] 



(2.53) 



(2.54) 



Unfortunately, we do not have a simple closed form for this integral. Note that the ratio 
P{s'^) = m^{^)/m\ is 0{e^). Using the one-loop potential ( p.26| ), one finds 

AQ T'^-' Att 



f{erPis'){l+0{q%s))} 



5.972 J Pis') (1 + 0(6)). 



(2.55) 



This is exponentially sensitive to 1/e because of the dependence on s. As with the correlation 



length, the overall numerical factor in (|2.55| ) is not useful unless s' is computed using the 
2-loop renormalization group. Phrased another way, the logarithm is 



In 



AQ 

,2rp3- 



-2 + e 



In(s') - Ine + 



In 



(ATiePit 



g2(s) 



7e - 



+ 0(6), 



(2.56) 
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and a 2-loop calculation of is needed to evaluate this result consistently through 0(e°). 
However, the need for the two-loop renormalization group is avoided if one computes the 
latent heat scaled by an appropriate power of the correlation length, 



4- Free energy difference, AF(T) 

We now consider the free energy difference between the symmetric and asymmetric states 
as the temperature moves away from Tp. It is convenient to use a dimensionless reduced 
temperature which we will define as 

(2.88) 

where Tq is the temperature below at which = is no longer metastable. In terms of 
the three-dimensional theory, this corresponds to using 

rri^^tml, (2.59) 

(neglecting higher order corrections in g\ and assuming that \t\ < 1). We shall not dwell on 
the computation of Tc — Tq, which involves details of the connection between the original 
four-dimcnsional and effective three-dimensional theories that are beyond the scope of our 
present study. 

To find AF(t) = T4ff(^, 0) — 14fr(i, (0)) (where (0) is the expectation in the asymmetric 
phase at temperature t), it is convenient to write 4> in terms of the rescaled field 

u = f{e)q4>/ii. (2.60) 



The one-loop potential then gives 



0.2302 rig T^-'s'-^D{t, e ^ 0) (1 + 0(e)) , (2.61) 
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D{t,e) = -[^^^) niin(2t^2_^^4-e^^^^4) ^ ^2.62) 



where 



4-e 

is the value of the right hand side at its non-zero local minimum. Note that D{0,e) = 1, 
D{1, e) = 0, and D{t, e) is smooth as e — > 0. 

Once again, the free energy density difference AF{t) is exponentially sensitive to e, but 
the free energy difference per unit correlation volume, 

CVmAi^W = I Dit, e ^ 0)(1 + 0(e)) , (2.63) 

is 0(e-2)_ 

5. Surface tension, a 

At the transition temperature T^, large domains of differing phases can coexist. Planar 
or nearly planar boundaries separating such domains will have a surface tension (or domain 
wall energy density) given, to leading order, by 







a = T d(j)J2V{(t)). (2.64) 



Using the change of variables (|2.60|) , this becomes 



X duj uj^Je - 2 cu2-^ + (2-e) uj^{l + 0{q{s))} (2.65) 
0.2987 ^ ( 1 + 0(e)) . (2.66) 



^. Bubble nucleation rate, Fn 

The bubble nucleation rate Fn is of order /i^exp(— 5e) where iSe is the action of the 
Euclidean bounce solution The rate increases from zero as the temperature drops 

below Tc. The relevant part of the action is the scalar contribution, 
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5e = 1 d^-'x + V{(j))) . (2.67) 

We shall look for 0(4— e)-symmetric bounce solutions. As discussed earlier, we may pa- 
rameterize the deviation of the temperature from Tc by writing = tm'^. The one-loop 
potential at A(s) = has the form V^^'^{(f)) = jj,'^"'^ v{q(j) / jj}'^^"^ ^t). The rescaling 

, l-e/2 

ct>^^ 0, x^—, (2.68) 

q ^iq 

moves the dependence on q out in front of the action, 

Se-^^J d'-'x (i(a0)2 + v{(l>, t)) . (2.69) 
q J 

We shall not numerically compute the extremum of the integral above, which depends on 
e and t. For small e and fixed t, however, the dependence of the action, together 

with fig. 1^, implies that the nucleation rate is larger than would have been found with an 
unimproved one-loop calculation. 

Instead of computing the action of an 0(4— e) symmetric bounce, one could also consider 
an 0(3) symmetric solution in analogy with the way we previously treated sphalerons. 

D. Testing the e-expansion when Ai <^ qf 

The motivation for exploring the e-expansion was the breakdown of perturbation theory 

in three dimensions. When Ai <^ qj, however, the unimproved loop expansion is a reliable 

method for studying the phase transition directly in three dimensions. This provides an 

opportunity to test the e-expansion in a regime where answers are already known by another 

technique. In this section, we shall assume that Ai ^ qf and expand our previous results 

in powers of xi = Xi/qf- To simplify the expansion, we shall also assume in this section 

that the initial three-dimensional couplings qf and Ai are small while Xi/qf is fixed. This 

corresponds to the formal limit ( [1.12| ) discussed in the introduction. When generalized to 

4— e dimensions, our limit becomes 

0^ A A 

— 0, with 4 < 1 and fixed. (2.70) 
e e qf 
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The couplings qf and Ai are still numbers of 0(e), but they are assumed to be small numbers 
times e. In particular, we shall ignore qf compared to eAi. 

In this limit, our results ( p.l8| ) and (|2.19| ) for q'^{s) and s"^ at A(s) = become 



q'is) ~ s^ql ~ ^ « 1, (2.71) 
aqf 

which can also be extracted directly from the x <^ 1 limit of the renormalization group 
equations ( p.4|) and ( p.lO| ). Since g^(s) is small, a one- loop calculation will be good to 
leading order in any dimension. 



Consider the one-loop result ( p.36| ) for the scalar correlation length in the asymmetric 
phase, where we had not yet assumed e to be small. In any dimension, this result will be 
valid to leading order in Xi/qf. Expanding in powers of e gives 

o s 1/2 / X l/e 



asym 



1 / 2 \ '-/'' / \ \ ^/"^ 

-1^1 I 4) [1.41748- 0.25701 e + 0.10441 e^- 0.03710 

+ 0.01760 - 0.00737e^ + O(e^) +0(xi)]. (2.72) 

We have chosen to factor out the leading power of e and the exponential dependence on 
l/e and then written the remainder as a series in e starting at e*^. Alternatively, we could 
have expanded In(^T) in powers of e as in ( [^.3(j| ). We shall generally present expansions in 



the form shown above, however, for little more reason than that we find the formulas more 
aesthetic. 

Consider the error which is introduced by truncating the series in ( |2.72| ) before setting e 



to 1. Truncating at 0(e°) and neglecting all terms proportional to positive powers of e yields 
an overestimate of C,c which is 14% larger than the correct three-dimensional value. If one 
instead keeps terms through 0(e), then the result is 6% too small. Clearly, the e-expansion 
for the correlation length is quite well-behaved, at least for small Xi/gf- For arbitrary Xi/gf, 
keeping terms to 0(e) requires a two-loop calculation of .^c and the three-loop renormalization 
group. Unfortunately, the three-loop /3-function for A is not yet known in either SU(2) or 
U(l) theory. 



Expansions like (|2.72| ) for some of the other quantities computed earlier are: 
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1 /g? V^V^^iV^^ 



T I eAi 



2.00462 - 0.86462 e + 0. 17588 - 0.09368 e 



eXiJ \aqi 



+ 0.03133 - 0.01698 + O(e^) + 0{xi)\ (2.73) 

-(2-e)/e 

1.36153 + 0.49372 e - 0.06630 - 0.00538 



cr = ngT 



\aqf 



- 0.00509 - 0.00021 + O(e^) + 0(xi)] (2.74) 

-(3-e)A 

0.14260 + 0.09226 e + 0.00692 - 0.00185 



AF(To) = ngT 



-(4-e)A 



- 0.00069 - 0.00028 + 0{e^) + 0{xi)] (2.75) 

0.23025 + 0.31089 e + 0.19796 + 0.10668 
+ 0.05882 + 0.03215 + O(e^) + 0(xi)j (2.76) 



The multiplicative errors made by keeping only terms through 0(e°) or O(e^) in the brackets 
are shown in table |. Note that the result for AF{Tq) is especially sensitive to the O(e^) 
term. 



observable ratio 


LO 


NLO 


asymmetric correlation length 


^asym 


1.14 


0.94 


symmetric correlation length 


^sym 


1.62 


0.92 


latent heat 


AQ 


0.77 


1.04 


surface tension 


a 


0.60 


0.98 


free energy difference 


AF(To) 


0.24 


0.56 



TABLE I. The ratio of the e-expansion resuhs, computing prefactors through leading order (LO) and 
next-to-leading order (NLO) in e, to the corresponding three-dimensional result when Ai <C (/i- 

For most quantities, we find that a calculation through 0{e^) gives agreement to within a 
factor of 2, and inclusion of O(e^) terms yields agreement within 10%. There is no guarantee, 
of course, that the convergence might not be worse when Xi/qf is large. The expansion of 



The series shown in Eqs. ( 2.72| - 2.76 ) are actually convergent at e = 1 because there are no singularities 



within |e| < 1 in the relevant formulas such as Eq. (2.35). This is only true at one-loop order; as shown 
later, higher order results do have singularities at e = ±1. 
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the free energy difference is notably worse than the other observables. In this particular case, 
the behavior of the series for the logarithm, ln(AF(To)/T^-'^), is rather different; next-to- 
leading order in e reproduces the correct three dimensional answer to within 17%. We have 
not shown the expansion of the baryon number violation rate because, due to the way we 
constructed its e-expansion, the result is trivially the same as the three-dimensional answer 
when Ai ^ qf. 



observable 


Co 


RGi 


Ci 


RG2 


C2 


RG3 


In^c, InAQ, InAF, In a 




0(6-1) 


0(6°) 


0(6°) 


0(6l) 


0(6l) 






0(6°) 


0(6°) 


0(6l) 


0(ei) 




InlB 


O(e-i) 


O(e-i) 


0(6°) 


0(6°) 


0(ei) 


0(ei) 


IuTn 




0(e-2) 


O(e-i) 


O(e-i) 


0(e°) 


0(e°) 



TABLE II. Calculations required for the e-expansion of various quantities discussed in the text. Co 
indicates a tree-level calculation, and C„ an n-loop calculation, at A(s) ~ 0; RG„ indicates n-loop renor- 
malization group flow. A computation through order e" requires all elements listed as e" in the table. 



As we have seen, various observables have different exponential dependence on 1/e and 
Consequently, determining their overall normalization can require computing to dif- 
ferent orders in perturbation theory. Table summarizes the loop order needed for both 
the renormalization group flow and the final calculation in order to compute the observables 
we have discussed to a given order in e. One generally requires using renormalization group 
evolution which is one higher order than the final calculation. (This reflects the fact that the 
renormalization group scale factor Ins is 0(l/e).) The only exceptions are ratios, such as 
those in the second line of table |I|, in which the leading dependence on s has been canceled. 



III. TWO-LOOP ANALYSIS 
A. Renormalization group flow 



The two-loop renormalization-group equations for the couplings have the form 
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s^ = eq'-/3oq'-P,q\ (3.1a) 
= e\- (aq^ + b q^\ + c A^) - (ko q^ + h g^A + ^2 g^A^ + A^) . (3.1b) 

as 

We shall later be computing the latent heat and so, as in the one-loop case, we also need 
the renormalization group equation for the mass, 

dlHj^ r 1 

s^ = - (/og^ + ^iA) + (nog^ + nig^A + n2A^) (3.2) 
as L J 

The new constants above have the following values: 



(4vr; 


Vi = 


-8N, 


SU(2) : {An] 


Vi = 


-§(544- 13iV) , 


(3.3a) 


{An 




-8 (45 + 7N) , 


{An 


fko = 


168 (73 - N) , 


(3.3b) 


{An 


fh = 


|(87 + 71Ar), 


{An 




-2 (455 - 71A^) , 


(3.3c) 


{An 


)'h = 


|(5 + 2iV), 


{An 


)'k2 = 


8 (5 + 2N) , 


(3.3d) 


{An 


th = 


-1(14 + 6iV), 


{An 


)'k,= 


-i(14 + 6iV), 


(3.3e) 


{An] 


^no = 


i(15 + 71Ar), 


{An] 


^rio = 


-527 + 71Ar, 


(3.3f) 


{An] 


^ni = 


f(l + iV), 


{An] 


^rii = 


16(1 + AT), 


(3.3g) 


{An] 


^n2 = 




{An] 


^^2 = 


-1(1 + 


(3.3h) 



The constants for the flow of A and have been extracted from the general results given 
in Ref. [^. As described in Appendix |B], the /5-function coefficients for the mass have been 



extracted from the results of Ref. [^, generalized to arbitrary A^ 



In contrast to the one-loop renormalization group flows, the two-loop trajectories are no 
longer independent of e. We will shortly derive the 0(e) correction to the one-loop solutions. 
Since we are eventually going to take e ^ 1, it is tempting instead to simply plug e = 1 
into the equations ( |3.1D and solve them numerically. This is not the same as solving for 
the expansion in e, truncating that expansion at next-to-leading order in e, and then taking 
e ^ 1. It is also a temptation that should be resisted. Fig. |lO| shows the result of solving 
the two-loop equations directly in three dimensions (e = 1) for the U(l) gauge theory with 
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= 1. The trajectories are qualitatively very different than the one-loop results of Fig. |^. 
This might at first seem like a disturbing problem for the e-expansion, but similar behavior 
occurs even in pure scalar theory, where the e-expansion is so successful. The line = of 
Fig. |10| corresponds to a complex pure scalar theory (or the xy model). The one- loop fixed 
point at A = 1/c has completely disappeared! The two- loop equations with e = 1 give no 
sign that the phase transition in the scalar theory is second order! This is equally true for 
the Ising model, which corresponds to N = 1/2, and whose critical exponents we reviewed 



in section |IB| . The lesson is that one must keep strictly to the philosophy of the e-expansion: 
take e to be small, compute results for physical quantities as a power series in e, truncate 
such series at some finite order, and only then send e to 1. 



g2 


► 




= ► 




> 


► 




> — 

■ *— 











Wr.- — ' X 




FIG. 10. Two-loop renormalization group flows in three dimensions (e ~ 1) for U(l) gauge theory with 
a single complex scalar. The right-hand figure shows an enlarged view of the region near the origin. 




FIG. 11. Two-loop renormalization group flows in three dimensions (e 1) for SU(2) gauge theory 
with a single scalar doublet. Flows below the heavy line run to A = +oo instead of curving back to A = 0. 
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For the sake of comparison, Fig. |TT] shows the analog of Fig. |T0| for SU(2) theory with 
a single Higgs doublet (A^ = 2). The violence to the behavior of the one-loop flows is 
less severe, but it is still qualitatively different. Of course, the three-dimensional behavior 
of the theory might really look something like Fig. |1^ instead of Fig. ^ with some flows 
running to A = +00; but the goal of our current study is to investigate the predictions of 
the e-expansion, in which case we must follow the philosophy outlined above. 

The 0(e) correction to q'^{s) from the two-loop RG equation is 



(3.4) 



where ql is the one-loop solution 



1 + 4(^^-1^ 
% 



and q^ is the O(e^) correction 



2/ \ /^l -e 4/ \ 



M J Wis) 



(3.5) 



(3.6) 



To find A(s), it is again convenient to solve the renormalization group equation for the ratio 
X = A/g^, 

dx 



ds 



-q^ [a+{b- (3o) x + cx^]-q^ [ko + {ki - (3i)x + k2 x^ + ks x^ 



(3.7) 



whose solution is 



xis) = Xa(s) + Xb{s) + O(e^), 



(3.8) 



where Xs,{s) is the one-loop result from ( p.ll| ) and 



^ ds' 
1 ~ 



■■'^{qtis') [ko + {k, - (3,)x^{s') + k,xl{s') + k,xl{s')] 
+ ql{s') [a + {h-^o)x.{s')+cxl{s')]]. 



(3.9) 



Here, 
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K{s) 



— qlis') [{f3o-b)-2cx^{s'] 

1 S' L 



2 In sec 



a 



2A 







2 In sec a 



with a and A given by Eqs. ( p.l3|) and ( p.l4|) , respectively. Solving x{s) 



s^ = sl + sl + 0{e' 



(3.10) 

for s now yields 
(3.11) 



where Sa is the one-loop result of ( |2.19| ) and 



^ (5a;a(sa)/9sa) aql{s^) 
Substituting this result into Eq. ( |H^ ) yields the value of g^(s) when A(s) vanishes, 

■ g^(^a)U0(e^) 



(3.12) 



e s 



4 dql{s^ 



(3.13) 



B. The two-loop potential 

The two-loop potential is one of the ingredients needed to compute physical quanti- 
ties, such as the correlation length, beyond the leading order in e. If the three-loop (3- 
functions were known, the two- loop potential could be used to obtain In^c to O(e^). Without 
three-loop /5-functions, the two-loop potential still allows one to determine ratios such as 
ln(A(5/^c '^) to O(e^) as summarized in Table 

The two-loop potential at A(s) = can be extracted in Landau gauge, with small mod- 
ifications, from the results of Ref. ||3^. In order to simplify the calculation, it is useful to 



note from Eqs. ( |2.32| ) and ( |2.33| ) that within the e-expansion the ratio of scalar to vector 
masses is small near the asymmetric ground state at Tc, 



2 2 

m m 



0{q') = 0{e). (3.14) 



M2 g202 

Hence, if one neglects corrections suppressed by additional powers of e, then the scalar mass 
may be set to zero when computing the two-loop potential. The details of extracting the 

47 



potential from Ref. [Q, as well as the result for general e, are given in Appendix Here, 
we shall simply quote the result when e — > 0, 



The coefficients Vi are: 

U(l): i;2 = i(9+iV), 

vi = -1(27 +4N) 



,6 J.4 



(47r) 



V2 In —7— + In ' 



;i + 0(e)). (3.15) 



r;o = ^[123+19A^ - {1-N)7i^ 



SU(2) : V2 = |(-61+A^), (3.16a) 

Vi = 286-4:N, (3.16b) 

I, vo = \[-19+19N - {2-N)7r^] (3.16c) 

- 297 



l + VSLi^j -y3^1n2 



TC 



where L{t) is Lobachevskiy's function, 

L{t) = — dx In cosx , 
Jo 



whose value at 7r/6 is 



L(7r/6) = 0.024617. 



(3.17) 



(3.18) 



C. The scalar correlation length 

As a sample computation using the two-loop potential, we shall discuss the scalar cor- 
relation length at the critical temperature. At one-loop order, the correlation length is 
completely determined by the curvature of the potential. 




(3.19) 



At higher orders, one must be more careful to examine truly physical quantities because the 
curvature of the potential is both gauge- and scheme- dependent. The physical correlation 
length ^ is determined by the pole position p = i/^oi the scalar propagator at zero frequency, 
which satisfies the dispersion relation 

p2 + + n(p) = , (3.20) 
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with po = 0. If one replaces Il{p) by 11(0), then the last two terms simply produce the 
curvature of the effective potential. The dispersion relation may therefore be rewritten as 

p' + nl + [Uip)-U{0)] = 0. (3.21) 

In the asymmetric state at T^, the curvature of the potential is order g'^fi^ while the mo- 
mentum dependent part of the self energy, Il{p) — 11(0), is order g'^p^fx/M ~ Locating 
the pole to leading order in g then requires only the one-loop effective potential, but at 
next-to-leading order one needs both the two-loop potential and a one-loop calculation of 

n(p)-n(o). 

At one loop order we computed the scalar correlation length in both the symmetric and 
asymmetric ground states. At two loops, the symmetric correlation length becomes prob- 
lematical in SU(2). To see this, consider the two- loop potential directly in three dimensions 
(using the e — 1 limit of the results of Appendix It takes the form 

f^V^'\^) - W ("^' + fzj + ^1 1' 1^ (^) + ^2 - 0{q') + 0{q') 0^ , (3.22) 

where Ci and C2 are constants. We shall focus on the 0^ term. The divergence as e — >■ 1 is 
the usual logarithmic ultraviolet mass divergence arising at two loops in three dimensions 
from graphs such as those illustrated in Fig. If we were doing perturbation theory 
directly in three dimensions, we would adopt a renormalization scheme that absorbed this 
divergence into the definition of a renormalized mass Tn^^^\ physical quantities would then 
have a finite relationship to m^^^ (rather than m?) as e — > 1. More troublesome than 



this ultraviolet divergence is the infrared behavior of the graphs in Fig. 12, which generate 



the 0^1nM ~ 0^1n(g0) term in the potential (|3.22|) . This term causes the curvature of 



the potential to diverge at = 0. For graphs (a) and (b) of Fig. [T^, this divergence is 
an artifact of the approximation where we neglected the scalar mass when evaluating the 
potential. For the non-abelian graph (c), however, it is unavoidable. As discussed in section 
P~F] , such infrared divergences are cut off by the non-perturbative physics in the symmetric 
phase responsible for the 0{q^n) mass gap for spatial gauge fields, and this scale will cut off 
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the logarithm in the two-loop correction to the scalar correlation length. This infrared cut- 
off of the (fP' ln(g0) term at scale where qcj) ~ g^/i also remedies a potentially embarrassing 



large > 19): = is then a local maximum instead of minimum, with a nearby local 
minimum at qcp ~ /iexp[— 0(l/g^)]. 



FIG. 12. Some two- loop contributions to the scalar mass that are logarithmically divergent in three 
dimensions. 

In the symmetric phase, the two-loop correlation length has a direct dependence on non- 
perturbative physics because this is the only source of an infrared cut-off. In the asymmetric 
phase, however, there is already a cut-off from the non-zero value of qcf). We shall therefore 
focus exclusively on the asymmetric phase, where the two-loop correction to the correlation 
length remains sensible even when e ^ 1. It would be interesting to compute multi-loop 
corrections to the symmetric phase correlation length within the e-expansion and investigate 
the extrapolation to e = 1, but we have not done so. 

To find the asymmetric phase correlation length, the first task is to compute the curvature 
of the two-loop potential at the asymmetric minimum 0c- Perturb around the one-loop 
solution by writing V = V^^^ + SV, 0c = 0c^'' + S(f)c, and ml = [m'^p]'^ + 5m? ^ and then 
linearize the equations V^(0c) = ^(0) and (9/90)V^(0c) = that determine the asymmetric 
state at the critical temperature. Solving for 0c and m^, one finds that 



feature of the perturbative result ( p.22| ) when ci is positive (which occurs for sufficiently 




(a) 



(b) 



(c) 




(3.23) 



(3.24) 



and the curvature of the potential at 0c differs from its one-loop value by 
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6Q 



asym 



-(3-e) 0, 



- d ^ f5V{^)-5V{0y 



(3.25) 



Plugging in the general results for 6V from Appendix |^ yields 



asym j'^^' ^ ^\ g7E 



1 + 



(47r)2 



Sie) + Oiq%s)) 



(3.26) 



where 



'S'u(i)(e) 



1 1 + f 1 



3 I g2 



(l-t)(l-t)' 



(iV_i)ei?(-e/2,e)- 
9 + AT 45 + TAT 



4 (5 - 4e + 
2 - e 



+ 



+ 



6e 



1 - 



(3.27) 



5; 



1 1 



SU(2)( 



2J \^ 4 



(-2 + A^ - (ll-4e) 3^^~''>/^ sin (e7r/2)) e5 (-e/2, e 
4 (20 - 35e + 21e'^ - Ae^ 



2 - e 



+ 9(ll-4e)(l-e)r(e) 
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_iV _ 7 (73 - A^) / 
6e V 



(3.28) 



Here, B{fi, v) is Euler's Beta function, and the function r(e) is defined in Eq. ( |A21| ) of the 
appendix. These results are indeed smooth as e ^ 1. The result for small e (which can be 
derived directly from the the e — ^ limit of Eq. ( ^.15D ) is 

1 



n„ 



5(6) = - |it;2 + |t;i-% [1 + 0(6)]. 



(3.29) 



FIG. 13. One- loop graphs contributing to n(p) — n(0). The black dot in the last graphs represents the 
one-loop counter term for wave-function renormalization. 



To complete the calculation of the correlation length, we need a one-loop calculation of 
n(p) — n(0) to order g^/x^. The graphs which contribute to n(p) — 11(0) are shown in fig. |T^. 
The calculation of these graphs may be simplified by remembering that m^/M^ ~ ~ 6 at 
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Tc. So and may be treated as small perturbations compared to when evaluating 
these graphs. In Landau gauge, one finds 



with 

= rig 



+ 



1 + 01 



-3n„ In 



M2 



1 + U 



2 gV ^ 

m? 



Plugging in the value of M 

6np 



using (|2.33|) gives 



S{e) 



1 



-1 



l-i + ^ 
6 12 



1 + 



2 ? P 



"2'^g 



l + 0(e,e^ 



(3.30) 

(3.31) 
(3.32) 

(3.33) 
(3.34) 



Inserting results ( p.26| ) and ( p.30|) into the dispersion relation ( p.21|) , the final result for the 
correlation length in the asymmetric phase is 

(3.35) 



asym 



D. Testing the e-expansion when Ai <^ qf 



As we did for the one-loop case, it is interesting to compare our two-loop corrections 
to the correct three-dimensional results in the perturbative regime Ai <^ qf. We first need 
the expansions of s"^ and q'^{s) in powers of Xi/qf at A(s) = 0. The expansion of s"^ can 
be obtained with some sweat by expanding the two-loop result ( p.8|) for x{s) and setting 
x{s) = 0. It is much easier, however, to return to the original two-loop RG equations ( p.4|) 
and solve them perturbatively in x = A/g^, recalling from the one- loop analysis (|2.71| ) that 



g^(s) ~ s'^qf ~ exi. One finds. 



q'{s) = s^qt-s^{s^-l)^Po + Oiexi), 

e 

2 4 
x{s) =x,- (s^-l) ^ (a + (6-/3o) x,) + {s^-lf^ ah - 

+ 0{x\ q\s) Ins). 



(3.36) 



(3.37) 
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Setting x{s) to zero and, as usual, assuming qf is small by taking the limit with xi 
fixed yields 

rW = -fl-^f''+-V0(x;)V (3.38) 



a \ 2a V a 
= !^ f 1 + ^ f2/?o - & - — U Oixl) ] . (3.39) 



aqf \^ ' 2a V ' " a 
Applying this expansion to our results ( p.27|) and (|3.28| ) for the two-loop curvature yields 



where Mo(e) = 1.41748 — 0.25701e + ■ ■ ■ is the previous series in Eq. (|2.72|) , and 

M^^^\t) = (0.23625 + 0.02625 A^) + (0.03591 + 0.02587 A^) 
+ (0.06521 + 0.03711 A^) e + (0.01481 + 0.02248 A^) 
+ (0.01249 + 0.02594 A^) + (0.01041 + 0.00964 A^) 
+ (-0.00035 + 0.01468 A^) + O(e^) , (3.41) 
^su(2) _ (_o.53374 + 0.00875 A^) e"^ + (-0.76655 + 0.00862 A^) 
+ (0.21722 + 0.01237 A^) e + (-0.26359 + 0.00749 A^) 
+ (0.21748 + 0.00865 A^) + (-0.22786 + 0.00321 A^) 
+ (0.2219 + 0.00489 A^)e^ + 0(e^). (3.42) 

The correct, three-dimensional result is 

M^^^\l) = 0.38119 + 0.17444A^, Mf^^^\l) = -1.24695 + 0.05815A^ . (3.43) 

The convergence of the e-expansion for the 0{xi) terms seems somewhat worse for the SU(2) 
case than the leading (xi)° terms that we examined previously.^ Nevertheless, in the SU(2) 



^^This can be traced to the pole of the Beta function in ( 3.28| ) at e = —1, which makes the radius of 
convergence equal to 1. This pole reflects the appearance of new, logarithmic two- loop divergences in five 
dimensions unrelated to those in four dimensions. 
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case with N = 2, when e — 1 the terms through 0(e°) in the coefficient Mi(e) of xi give a 
result differing from the three-dimensional result by 12%, and the sum of terms through O(e^) 
err by only —9%. After that, higher-order terms give larger and larger contributions, which 
is the signal to stop. So, either the 0{e^) or O(e^) results give reasonable approximations to 
the correct answer. Keep in mind that, since we have assumed Xi is small, we are currently 
testing the e-expansion on a small correction to the total result. 



E. The latent heat as a test for Xi/qf > 1. 

We shall now test the e-expansion in the range Xi/qf ^ 1, where unimproved pertur- 
bation theory fails, by computing a sample quantity to next-to-leading order and checking 
whether the correction to the leading-order result is large or small. Recall from the discus- 
sion in section [II C| that, without three-loop RG equations, we cannot consistently compute 



next-to-leading order corrections to the prefactors of quantities such as the correlation length 
or the latent heat. But we can compute the correction to the ratios such as ^I^^^AQ with 
a purely two-loop calculation using the two-loop renormalization group. In this section, we 
do a next-to-leading order calculation of ^"^^^AQ to test the behavior of the e-expansion. 

We first need to compute AQ by using (p.52|) and linearizing about the one- loop result. 
We won't derive the two-loop result for arbitrary dimension but shall instead simply use the 
e — *• limit ( p.l5|) for the two-loop potential. One finds 



Z/2 S^-' q\s) ' ' \ (477)2 Tig V27 - 9 

and ^iy^AQ follows from the results (ICTI) , ([Q^ ) and (lOSi) for ^ 



asym) 



e^y^AQ 4! P{s^) q\s) 



0{^) . (3.45) 



Here, P{s'^) is the two loop version of the ratio ( p.54| ) of the running mass at the initial and 
final scales. 
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exp< — 



e Ji s" 



+ no g^(s') + rii g^(s')A(s') + n2 A^(s') 



= Paisl) [l + S^ + Oie% 
where Pa is the one-loop mass ratio, 

Paisl) ^ expj-- qlis') [/o + h x.{s')]j 

and 6i is the two-loop correction, 



(3.46) 



(3.47) 



Si = (3a) - - r ^ ^o g^(s') + h [ql{s')x^{s') + x.{s')ql{s')] 



(3.48) 



Finally, we need the expansion ( p.l3| ) of g^(s) at A(s) = 0, or 

1 1 



where 



(e - q 



+ 



gb(ga) 

?a(^a) 



(3.49) 



(3.50) 



Inserting these expansions of g^(s) and Pis'"), plus that of the scale s% in Eq. (|3.11| ), into 
the result (|3.45| ) yields 



Z/2 



2a qtis 



[l + 5tot + 0(e2)] 



where 



(47r) 



(3.51) 



(3.52) 



The relative size of this 0(e) correction is shown in Fig. [1^ for the minimal standard 
model (in our approximation of ignoring the Weinberg mixing angle). The correction varies 
between roughly ±30% for (zero-temperature) Higgs masses up to 150 GeV. This suggests 
that the e-expansion is tolerably well-behaved for these masses. For larger masses, the 
correction does not increase indefinitely but is bounded by 80%. This gives hope that the 
e-expansion may be a useful qualitative description of the phase transition even for these 
larger masses where it does not work as well quantitatively. 
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FIG. 14. The relative size of the next-to-leading order correction to ^^^y^-^AQ in the e-expansion. The 
values are given as a function of the (tree-level) zero-temperature Higgs mass in minimal SU(2) theory 
{N = 2) with g = 0.63. 

The comparatively small size of 6tot becomes more impressive when one examines the 
size of the four separate terms in (|3.52|) . At mh(0) = 80 and 250 GeV, they are 

5tot( 80 GeV) = -0.45 + 4.88 - 0.94 - 3.79 = -0.30 , (3.53) 
5tot(250 GeV) = -4.42 + 34.83 - 20.00 - 9.66 = 0.75 . (3.54) 

These large cancellations clearly underscore the importance of examining physical quantities, 
rather than unphysical ones such as g^(s) or s% when testing the e-expansion. 

IV. MANY SCALAR FIELDS 

In this section, we examine two calculations related to theories where the number of 
complex scalar fields, N, is large. As we shall see, the e-expansion is less well behaved in 
such theories. 

A. Critical A^c 

As noted in the introduction, a non-trivial infrared stable fixed point exists in both U(l) 
and SU(2) theories if the number of scalar fields is sufficiently large, N > N^. As e — > 0, 
there can never be an infrared-stable fixed point if the gauge coupling is asymptotically 
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free, since then the one-loop renormahzation group equation (|2.4a|) has no non-trivial zero. 
A critical number A*",, of scalar fields exists for SU(2) only because jSo becomes positive for 
> 88, thereby permitting a non-trivial fixed-point for sufficiently large A^. If there is a 
fixed point, it must have 



Q' = i + Oie'). (4.1) 

A fixed point must also have dx/ds = 0, where this derivative can be obtained from Eq. ( p.Tj) . 
As A^ approaches Ac from above, the stable fixed point A and the tricritical point B shown in 
Fig. ^ merge. This provides a simple way to compute A^^, for when the two fixed points touch 
not only must the right-hand side of the RG equation ( p.7| ) vanish, so must its derivative 
with respect to x. In particular, the combination (1 — \xd/dx) must give zero when applied 
to the right-hand side, which implies that 

(a + i(6 - (3o) X,) + ql {ko + Pi)x, - \h a^') = . (4.2) 

This has the solution 

Po - Po-o \ Po - {/Jo-byj 

Plugging back into the RG equation dx/ds = gives 

Arc = ArW+eAr«+0(e2) (4.4) 

where A^^^^ is the value of A^ for which A = 4ac — (h — I3qY is zero and A"^^^) is 



Inserting the coefficients from ( ^Tsf) yields the results 



N=N, 



(0) 



(4.5) 



A^c^^i) (4-e) = 182.95 - 242.67 e + 0{e^) , (4.6) 
A'SU(2)(4_e) = 718 -990.83 e + 0(e2). (4.7) 

As noted in the introduction, these series clearly can not be used to obtain a reliable estimate 
of A'c when e = 1. 
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B. Tricritical slope X/q^ 



The goal of this section is to compute, for large A^, the ratio \/q^ of the tricritical point 
separating the domains of first and second order phase transitions, as shown in Fig. ^ The 
lowest-order term in the e-expansion is given by the N ^ oo limits of the expression (p.l7|) 
for X- . The first order correction is easily determined by solving for the zero of the two loop 
renormalization group equation ( |3.7| ) to order e. One finds 

A 



y2 



54 - 126e + 0(e^)] ^ + O ) . (4.8) 



This e-expansion is also poorly behaved. We shall see below that, though the qualitative 
dependence on is correct, the leading-0(e'') coefficient of 54 differs from the correct three- 
dimensional coefficient by a factor of roughly five! It is important to emphasize, however, 
that the e-expansion alerts one to its own failure by producing a next-to-leading order 
correction that is significantly larger than the leading-order result when e = 1. 

Even the qualitative result that the critical A/g^ scales like is at odds with a large 
analysis in three dimensions by Jain and Papadopoulos ||25|. One of the main purposes 
of this section is to redo their analysis and show that the e-expansion result is qualitatively 
correct. We shall closely follow their analysis, though our final result is different. It should 
be emphasized that large in the present context means a large number of scalars, not the 
replacement of the gauge group by SU(A^). 

Before focusing on three dimensions, we note a qualitative feature of the e-expansion 
result which will help guide us in selecting the orders of parameters of interest in the three 



dimensional case. The tricritical point determined by ([4.8|) approaches the A = axis as 
A^ —>■ oo. Recalling that the value of at the tricritical point scales like 1/A^, the prediction 
for A at the tricritical point is 0(1/A^^). 

Now consider the three dimensional theory (with either U(l) or SU(2) gauge group). 
Taking a hint from the e-expansion, examine the regime where is order 1/A^ and A is 
order We want to compute the effective potential in the large- A^ limit and determine 
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whether its variation with yields a first- or second-order phase transition. Since the 
effective potential is a function of 0, evaluating its large- limit depends on how (f) scales 
with A^. To get a hint, recall the results of naive one-loop perturbation theory. This produces 
a potential of the form 

y (0) ~ im202 _ g3^3/2^3 + ^ A/i 0^ (4.9) 

where the constant k is 0(1). In the asymmetric phase at the transition, all three terms are 
of the same order, which implies that 



, O(VAr), (4.10a) 
A 

^ X^icj)'^ = 0{l/N) . (4.10b) 

In this range of the vector mass M is large compared to m, 

= = 0(A^°) . (4.11) 

Using these parameters, we may easily compute the three-dimensional effective potential 
to leading order in A^. The dominant graphs are shown schematically in Fig. ^ and are 
vector loops decorated by scalar rings. Any other type of graph is sub-leading in A^. Hence, 
in Landau gauge, the leading-order potential (neglecting 0-independent terms) is simply 

V^(0) = imV + ¥V0' + %/ ^ Hp' + M^ + UM). (4.12) 

Here, Hv(p) is the scalar contribution to the one- loop vector self-energy. The dominant 
momenta flowing through the ring graphs are of order M. Consequently, in the self energy 
Hv(p) the scalar mass m? is small compared to and may be neglected. The leading 
contribution to Hv(p) is then easily computed to be 

H.(p) = ^ArgV- (4.13) 
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FIG. 15. Dominant graphs in large N for X/q^ ~ 



We now need to evaluate the integral in ( [4.12|) . This integral must first be regulated and 
renormalized. Two subtractions are needed to render the integral finite, one 0-independent 
and one proportional to 0^. The resulting finite, renormalized integral is 



where 



(2^ 



In 



^ p2 ^ nv(p) , 



M2 



p2 + Il^{p) 



-\M'f{z) 



(4.14) 



-3/2 



(2-4)Vi^S 



ec \/z 



|.--2)l„(|l+.- 



(4.15) 



and 



_ 1024 _ 1024 02 
The resulting potential can then be written in the form 



(4.16) 



lm202-|n,0V(lO240^) + iT^0' 



(4.17) 



where = M/{Nq^ij) = cj)/ [NqjjL^/'^). The condition for a first-order phase transition is 
that there exist a mass fh and field 7^ such that V{(f)) = V'{(j)) = 0. By applying 
(1 — 1(09/90)) to Eq. ( |4.17|) , one may derive the constraint 



d 



(0/(102402) 



1 NX ~ 



(4.18) 



The left-hand side is plotted in Fig. 0, and the right-hand side is depicted by the dashed 
line. Because the left-hand side is concave downward, it is clear that the maximum value of 
N\/q^ which satisfies the constraint is 
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An, 



0/(102402 



0=0 



96 n„ 



(4.19) 



This is precisely the condition that the coefficient of cj)^ vanish in the expansion of V{(f)) 
about = 0, which is the standard mean-field criterion for a tricritical point. 



l/27t -- 



NA.$/4ngq2 




d((()f)/d$ 



FIG. 16. Graphical representation of solution to Eq. (4.18) 



Before comparing this result to the e-expansion result of ( |4.8| ), one must convert the 
definitions of A and g^, which are scheme-dependent. In our three-dimensional analysis, A 
and have simply referred to the (dimensionless) bare couplings Aq and gg, which do not 
require any ultraviolet renormalization. In the e-expansion, however, our A and were 
defined by minimal subtraction in four dimensions, so that 

1 



Ao = A + -(ag^ + hq'\ + cY) + O(g^), 



ql = q^ + -f3oq' + 0{q' 



(4.20) 
(4.21) 



Amusingly, at leading order in 1/A^, (Ao,go) ^^e tricritical point is twice (A, g^) and so 
the ratio A/g^ is not changed by the redefinition. 

The three-dimensional result ( [4.19[ ) has the same N and Ug dependence as does the e- 
expansion result (|4.8| ), but the coefficient differs from the ffist term of the e-expansion by a 
factor of roughly five. As discussed before, the ill-behaved e-expansion ( [4.8| ) clearly cannot 
be used to obtain a meaningful estimate of the overall coefficient at e = 1. 

The preceding discussion of the scheme-dependence of A and g^ brings up an important 
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point. The ratio A/g^ at the tricritical point is not a physical quantity; it is scheme- 
dependent and will change under generic finite renormalizations. At the very end of section 
p.11 E| we discussed the dangers of testing the e-expansion with unphysical quantities. One 
should really test the e-expansion in the case at hand with something physical. A good choice 
would to consider the RG fiow line that connects the Gaussian fixed point at A = = 
with the tricritical point. The slope A/g^ of the line at the Gaussian fixed point is invariant 
under finite, perturbative redefinitions of the couplings because the couphngs are arbitrarily 
small near the Gaussian fixed point. Within the e-expansion, the leading-order result for 
this slope is no different from the result 5Ang/N at ( [4.8|) for the tricritical point, and so 
again differs from the three-dimensional result by a factor of roughly five. We have not 
computed the expansion of this physical ratio A/g^ to higher orders in e. 

Contrary to our results here, Jain and Papadopoulos saw no evidence of a first-order 
phase transition for any value of A/g^ of order They missed it because they did not 

look for the transition in the proper range of and m^, as given by Eq. ( [4.10a] ). 

V. CONCLUSIONS 

We have applied e-expansion techniques to the study of first-order electroweak phase 
transitions. The e-expansion is not always reliable, but it provides its own diagnostic for 
failure. A comparison of whether sub-leading terms in the e-expansion are larger or smaller 
than the leading term provides an a posteriori indication of whether the e-expansion is likely 
to be a reasonable approximation for the problem at hand. 

The e-expansion has had some spectacular successes, such as the computation of crit- 
ical exponents in the Ising model. Assessing its utility in gauged Higgs theories is more 
complicated. By comparing with large calculations, we have found that the e-expansion 
works qualitatively but does not provide quantitatively useful predictions at e = 1 when the 
number of scalar fields is large. Nevertheless, in this case, we reassuringly also found that 
the sub-leading terms in the e-expansion were large compared to the leading one; so the 
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expansion diagnosed its own failure. 

In contrast, our results for real problems of interest at small are much more encour- 
aging. Several tests of the e-expansion in the perturbative regime X <^ g^, corresponding to 
very light Higgs masses, showed that the expansion is quite well-behaved. However, what 
we really want to know is how well the e-expansion does when A > (7^ and conventional 
perturbation theory fails. We have made one non-trivial test in this case by computing to 
next-to-leading order in e the latent heat of the transition (scaled by the correlation length 
squared). For Higgs masses below 100-200 GeV, this e-expansion does not seem too bad. In 
contrast, we know that for some quantities conventional perturbation theory fails badly for 
Higgs masses as small as 35 GeV. We have only made one such test of the e-expansion for 
A > (7^, however, and more are certainly called for. It's always possible that we happened 
to pick a particular quantity that is better behaved than most. 

If, based on the above evidence, we tentatively assume that the e-expansion is a reason- 
able approximation, then we can glean insight from a number of relatively simple leading- 
order calculations. In particular, leading-order calculations predict that the transition is 
stronger than that given by conventional one-loop (ring-improved) perturbation theory when 
the Higgs mass is below roughly 130 GeV. However, in contrast to intuition based on the 
tree-level relation between the sphaleron mass and the Higgs expectation (0), the rate of 
baryon number non-conservation in the asymmetric phase is always larger. Even if, in the 
end, one finds that these calculations do not extrapolate well from 4— e to three dimensions, 
there is still an important lesson to extract. If the equilibrium properties of the electroweak 



phase transition are studied on a lattice [35|, and if the phase transition is found to be 



considerably stronger than expected, this does not necessarily make electroweak baryogen- 
esis more viable in the model studied; an unexpectedly strong transition is not necessarily 
one with an correspondingly small rate of baryon nonconservation at the completion of that 
transition. 

Based on the success of our main test of the e-expansion, it is possible that the expansion 
will provide a useful quantitative as well as qualitative tool for cases of interest to electroweak 
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baryogenesis. The particular application to the baryon nonconservation rate, however, forces 
one into awkward technical contortions related to the intrinsically three-dimensional nature 
of the sphaleron. More study will be needed to see if the suggestions we have made will 
really produce a successful method for computing this rate to next-to-leading order in the 
e-expansion. 

This work was supported by the U.S. Department of Energy, grant DE-FG06-91ER40614. 
We thank Sergei Khlebnikov and Larry Sorensen for useful discussions. 
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APPENDIX A: THE TWO-LOOP POTENTIAL 



The two-loop potential can be extracted from Ref. if one specializes to U(l) or SU(2), 
and changes the overall group factor for each diagram to generalize to an arbitrary number 
of complex scalar degrees of freedom. Setting A to zero and ignoring compared to 



M^, the relevant terms in eqs. (5.4-5) of Ref. become 



(Al) 



where 



Vsy = ngg2{|iVA(0, 0,z)-z Biz, z,0)} , 
V^^u(2) _ q'^{_2A{z, z, z) + 6S(^, z) - 3A(0, 0, z)} . 



(A2) 
(A3) 
(A4) 



The functions A, 5, A and S are given in eqs. (5.8-15) of Ref. and are: 



A(0,0,z) = -^/(z,0,0)-^J(z), 

B{z, z, 0) = 3/(^, z, 0) - ^/(z, 0, 0) - 2 el{z, z, 0) + ^^(^, z) + 3 J(z) 

99 - 3 - 

A(z, z, z) = ——zl{z, z, z) + -^zl{z, 0, 0) + 18 zel{z, z, z) 

+ 9 J{z, z)~^- + 4) -eJ{z, z) - 61 zJ{z) , 

S(z, z) = ^J{z, z) + [fc^ - ^) J(z, z) - 9 zJ{z\ 



(A5) 
(A6) 

(A7) 
(A8) 



where 



z = q 



212 



(A9) 



and we have defined J{x) in terms of their J{x) by^ 



J{x) = J{x)/{47if 



(AlO) 



^^There are typographic or notational errors in Ref. 1^ ; equations (5.8-15) are correct only if each exphcit 
factor of J{w) in those equations (but not J{v,w) or J{v,w)) is replaced by J{w). Also, the factors of /i^*^ 
in their Eqs. (3.4-5) should be eliminated. 
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Their definition of e, written as e in this appendix, is related to ours by 



e = e/2 . 



(All) 



The functions J, J, /, and / are defined in Eqs. (3.1-5) of Ref. The special cases of 
these functions which are needed may be extracted from eqs. (4.13), (4.16-7) and (4.19) of 



Ref. and are: 

J{z, z 
1(^,0,0 
J(z,2,0 

J(z, z, z 

J{z, z 
1(^,0,0 
J(z,2,0 
/(z, z, z 



(47r)4 \^47r/i2 ) 

^2 / ^ \ 



r(-i+e) 



(47r)4 \ 47r//2 



^(-l+e)^ 

■r(e)r(-i+2e)r(i-e) 



(47r)4 Utt/x^ 



-2e 



r(e)r(-i+e) 



(47r)4 \^47r/i2^ 



1 -2e 



-2e 



(47r)4 \^47r/i2 ) 
J{z,z) + 2e~hJ{z) , 
l{z,0,0)-e-'j{z), 
l{z,z,0) -2e-Kj{z) , 
I{z, z, z) — ?)e~^J{z) . 



I r(2e) r(e)r(-l+e) + /(Vs^, 0, 0) sin(e7r) 



(A12) 
(A13) 
(A14) 
(A15) 

(A16) 

(A17) 
(A18) 
(A19) 
(A20) 



We have also used the limiting values J(0) = J{z, 0) = J{z, 0) = C{z,Q) =0 to simplify the 
results of Ref. ||3^. The function r appearing in I{z, z, z) is 

-•i At /4\^" / 1 3 r 



r(2e) 



-)2e 



(t2 + 3)^ 



^^^'^'2' 2' 3 



and has the limits 



r{x) = 1 + X 



1 + 



(l + ln2-iln3) 



^(1) = 2Sinh-^(l/V3). 



(A21) 



+ 0{x^) , (A22) 
(A23) 



Here, L{t) is Lobachevskiy's function, defined in Eq. ( ^.171 ). Note that to get a result for 
general e, it is important to use formulas (3.3) and (4.19) of Ref. p3 where the limit e ^ 
has not yet been taken. 
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Putting everything together, and switching from e to our canonical definition of e, gives 
the following results for U(l) and SU(2), 

3-e 



.yd) 



(47r) 



1 -e 
1- 



r(-i+ie)r(ie)-ir(-i+ie^^ 



+ 



2-e 



r(i-ie)r(|e)r(-i+6 



2-e 



+ {9 + N) {l-l)n-lHe) 



q 



4:71 jj"^ ^ 

2 72 \ 2-e/2 



+ 



9 + N 45 + 7iVl /gV^ 



6e 



(A24) 



A'' '^SU(2) — P '^SU(2) 



(4vr) 



+ ^ (11 - 4e) r(e)) r(-l+le) r(ie) 
+ 3(3^^)/27rr(-l+e) + (7 - 8e + 2e^) Ti-l+^ef 



+ 



2-e 
2- AT 
2-e 



r(i-ie)r(ie)r(-i+e) 



2-e 



+ (-6i + iv) (7-^)r(-i+| 



9 



2 72 \ 2-e/2 



+ 



-61 + iV 7(73-A^) 



6e 



47r/i2^ 

2 j;2 \ 2-v 



47r/x2 



X 



(A25) 



APPENDIX B: THE TWO-LOOP SCALAR MASS /^-FUNCTION 



Ref. gives standard model results for the /3-function for the scalar mass. One may 
easily generalize these results to the case of an arbitrary number of scalar fields, just as was 
done for the two-loop potential in Appendix ^ The only change necessary is to insert the 
appropriate dependence into the group factor of each graph, which changes Eqs. (5.2) 
and (5.4) of Ref. Q to 



12 



I{H, H, H) + l(2iV-l) I{H, G, G) 



+ 



A r 



J{H, H) + |(2Ar-l) J{H, G) + \{AN^-l) J{G, G) 



(Bl) 
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sv 



A{H, G, Z) + (iV-1) A{G, G, Z) - 2Z B{Z, Z, H) 
+ C{Z,H) + {2N-1) C{Z,G)\, 



where 



2 D 



and Z = q 



2 J 2 



(B2) 



(B3) 



Following the derivation of Ref. and inserting the appropriate Landau gauge values of 



the scalar anomalous dimension 33 



d(j) 
^ dii 



■^(l) + ^(2) + o(A^g^) 



(47r)2 7« = -3ngg% 

\4^,(2) _ 1 I 1^ \2 , 1 



(B4) 

ajji 
with 

(B5) 
(B6) 
(B7) 

then yields the /3-functions given by our (|3.3|) . Because we are interested in non-zero values 
of A, evaluating the resulting expressions is slightly more complicated than was the case for 
Appendix 0; however it is also easier since one may restrict attention to the e — > limit. 



(4vrr 7|;(^) = (iV+1) + I (lliV+9) g^ 
(4^)' ifim = Ts (N+l) X' + I (lliV-533) q' 
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